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HYMO: 

PROBLEM-ORIENTED  COMPUTER  LANGUAGE 
FOR  HYDROLOGIC  MODELING 

Users  Manual 

By  Jimmy  R.  Williams1  and  Roy  W.  Hann,  Jr.2 


HYMO  (7)3  is  a  problem-oriented  computer 
language  for  modeling  surface  runoff  and  sedi- 
ment yield  from  watersheds.  The  language  is 
called  HYMO  from  the  words  "hydrologic 
model."  HYMO  was  designed  for  planning  flood 
prevention  projects,  forecasting  floods,  and  re- 
search studies.  It  consists  of  a  main  program  and 
16  subroutines  written  in  FORTRAN  IV,  but  it 
can  be  used  by  hydrologists  with  little  knowl- 
edge of  computer  programing.  The  language 
provides  17  commands  for  the  hydrologist  to 
use  in  any  sequence  for  application  to  any 
watershed. 

HYMO  was  designed  to  transform  rainfall 
data  into  runoff  hydrographs  and  to  route  these 
hydrographs  through  streams  and  valleys  or 
reservoirs.  It  will  also  compute  the  amount  of 
sediment  produced  by  a  storm  at  any  point  on  a 
watershed.  It  will  be  useful  to  research  hydrolo- 
gists in  studying  the  effects  of  watershed  and 
storm  characteristics  on  the  flood  hydrograph. 
HYMO  is  also  a  good  research  tool  for  testing 


hydrologic  procedures;  for  example,  a  new 
flood-routing  method  could  be  added  to  HYMO 
and  tested  easily,  because  the  inflow  hydro- 
graphs  and  the  rating  curves  are  available  in  a 
HYMO  program. 

HYMO  is  flexible.  Present  hydrologic  proce- 
dures can  be  modified  or  deleted,  and  other 
hydrologic  procedures  can  be  added  by  hydrolo- 
gists familiar  with  FORTRAN  IV  programing. 
Adding  a  new  command  simply  requires  the 
addition  of  a  new  subroutine. 

HYMO  is  efficient,  practical,  and  generally 
applicable.  HYMO  programs  can  be  written  and 
the  results  interpreted  by  hydrologists  who  have 
no  conventional  computer  programing  experi- 
ence. The  hydrologic  procedures  used  in  HYMO 
are  practical  —  required  inputs  are  easy  to  obtain 
for  most  watersheds. 

HYMO  was  written  for  the  IBM  360-65 
computer,  but  it  could  be  run  on  an  IBM  1130 
with  little  modification.  The  storage  require- 
ment is  about  73  K. 


OPERATION  OF  HYMO 


HYMO  consists  of  a  main  program  and  16 
subroutines.  The  HYMO  card  deck  is  set  up  in 
the  following  order: 


'Hydraulic  engineer,  College  Station  Area,  Southern  Region, 
Agricultural  Research  Service,  U.  S.  Department  of  Agri- 
culture, Riesel,  Tex. 

"Head,  Environmental  Engineering  Division,  Civil  Engineer- 
ing, Texas  A&M  University,  College  Station,  Tex. 

italicized  number  in  parentheses  refer  to  items  in  "Liter- 
ature Cited"  preceding  the  appendix. 


1.  Main  program. 

2.  Subroutines. 

3.  A  data  card  containing  the  number  of 
commands  in  the  command  table. 

4.  A  data  card  containing  the  ZALFA  array. 

5.  Seventeen  data  cards  containing  the  com- 
mand table. 

6.  The    users   program   deck   consisting  of 
program  and  data  cards. 

A  printout  of  the  main  program,  subroutines, 
ZALFA  array,  and  command  table  is  given  in 
the  appendix. 
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The  main  program  reads  the  command  table 
and  then  calls  the  HONDO  subroutine  (2)  to 
read  a  program  data  card.  Subroutine  HONDO 
determines  the  command  name  and  number  by 
comparing  columns  1  through  20  of  the  program 
data  card  with  the  command  table.  Then 
HONDO  determines  individual  data  items  by 
comparing  columns  21  through  80  of  the  pro- 


gram data  card  with  the  Z ALFA  array.  The  data 
are  placed  in  an  array  and  returned  to  the  main 
program.  Based  on  the  command  number,  the 
main  program  calls  the  proper  subroutine  to  do 
the  desired  calculations.  When  the  calculations 
are  complete,  control  is  returned  to  the  main 
program,  and  HONDO  is  called  again  to  read  the 
next  program  card. 


HYDROLOGIC  PROCEDURES  USED  IN  HYMO 


Trie  procedures  used  in  HYMO  were  selected 
because  of  their  accuracy,  general  applicability, 
practicality  of  inputs,  and  computational  effi- 
ciency. For  most  watersheds  the  input  is  easy  to 
obtain,  and  the  procedures  produce  reasonably 
accurate  results  without  excessive  computer 
time. 

Hydrograph  Computation 

When  flood  routing  is  performed,  a  watershed 
is  divided  into  many  small  areas  according  to  its 
hydraulic  characteristics.  The  hydrographs  from 
these  areas  must  be  estimated,  since  streamflow 
measurements  are  seldom  available.  A  procedure 
for  computing  unit  hydrographs  was  developed 
previously  (4).  A  modification  of  this  procedure 
is  used  in  HYMO.  Unit  hydrographs  are  divided 
into  three  parts  for  computation  (fig.  1).  From 
the  beginning  of  rise  to  the  inflection  point,  t0 , 
the  hydrograph  is  computed  by  the  two- 
parameter  gamma  distribution  equation 


Qp 


(n-l)  (i-w)(f/fp-i) 
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(1) 


Figure  1.  Dimensionless  unit  hydrograph. 


where   q  =  flow  rate  in  cubic  feet  per  second  at  time  t, 

qp  =  peak  flow  rate  in  cubic  feet  per  second, 

tp  =  time  to  peak  in  hours, 

and       n  =  dimensionless  parameter. 

From  t0  to  t1  (tt  =  t0  +  2K)  the  hydrograph 
is  computed  by  the  recession  depletion  equation 
tn-t 


K 
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(2) 


where  q0  =  flow  rate  at  the  inflection  point, 


0  =  time  at  the  inflection  point, 
and      K  =  recession  constant  in  hours. 

From  ty  to  oo  the  recession  depletion  equa- 
tion becomes 


q  =  qxe 

where  ql  =  flow  rate  at  tY , 

and    Ki  =  3K  =  second  recession  constant. 


(3) 


The  dimensionless  shape  parameter,  n,  is  a 
function  of  K/tp,  as  shown  in  figure  2.  The  peak 
flow  rate  is  computed  by  the  equation 


BAQ 


Qp 


(4) 


where  B  =  a  watershed  parameter,  a  function  of  n  as 
shown  in  figure  3, 
A  =  watershed  area  in  square  miles, 
and      Q  =  volume  of  runoff  in  inches. 

Therefore,  the  entire  unit  hydrograph  can  be 
computed  if  K  and  tp  are  known.  K  and  tp 
can  •  be  determined  by  hydrograph  analysis  (4) 
for  gaged  watersheds.  To  compute  K  and  tp  for 
ungaged  watersheds,  HYMO  uses  the  equations 


K  =  27.0A0.231SLP-0.m(L/W)0.l24 


(5) 
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and        tp  =  4.63^0.4225^-0.46^/^0.133  (6) 

where  SLP  =  difference  in  elevation  in  feet,  divided  by 
flood-plain  distance  in  miles,  between 
watershed  outlet  and  most  distant  point 
on  the  watershed, 

and    L/W  =  watershed  length-width  ratio. 

Storm  hydrographs  are  computed  by  con- 
volving unit  hydrographs  with  incremental 
source  runoff.  To  compute  incremental  source 
runoff,  the  mass  rainfall  curve  is  broken  into 
equal  time  increments,  and  the  Soil  Conserva- 
tion Service  (SCS)  rainfall-runoff  relationship 
(5)  is  applied.  The  SCS  rainfall-runoff  relation- 
ship is  expressed  in  a  set  of  numbered  curves. 
The  SCS  National  Engineering  Handbook  (3) 
provides  detailed  instructions  for  selecting  the 
proper  curve  number. 

Hydrographs  computed  by  this  procedure 
compared  closely  with  measured  hydrographs 
from  34  watersheds  located  in  Texas,  Oklahoma, 
Arkansas,  Louisiana,  Mississippi,  and  Tennessee. 
The  watershed  areas  ranged  from  0.5  to  25 
square  miles. 

Flood  Routing 
Streams  and  valleys 

The  variable  storage  coefficient  (VSC)  flood- 
routing  method  (5)  was  selected  for  HYMO. 


Figure  2.  Relationship  between  dimensionless  shape 
parameter  and  recession  constant/time  to  peak. 
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Figure  3.  Relationship  between  dimensionless  shape 
parameter  n  and  watershed  parameter  B. 


The  VSC  method  has  been  revised  (6)  to 
account  for  the  variation  in  water  surface  slope 
during  a  flood.  The  revised  VSC  method  is  about 
as  accurate  as  the  implicit  method  (7)  and  has 
the  general  applicability  of  simpler  storage 
methods.  Although  an  iterative  solution  is  used, 
the  VSC  method  requires  little  computer  time 
and  is  free  of  convergence  problems. 
The  VSC  routing  equations  are 


O, 


=  C2[la  +(^--iJ  o}. 


c,  = 


c,  = 


2  A  t 

2  A  t 
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1800^  +  VQl) 
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L  x  SLP0 
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(10) 


T?  = 


1800 /T!,  +  V0J 


X 


L  x  SLP0 


LxSLP0+Du  D0„ 


(11) 
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In  these  equations  subscripts  1  and  2  refer  to 
the  beginning  and  end  of  the  time  interval  a*;  the 
units  are  cubic  feet  per  second  for  flow,  hours 
for  time,  feet  per  second  for  velocity,  and  feet 
for  length  and  depth.  The  symbols  are  defined  as 
follows: 


r 
1 

=  inflow  rate. 

(J 

=  outflow  rate. 

h 

2      =  average  inflow  rate. 

C 

=  storage  coefficient. 

T 

=  travel  time  through  the  reach. 

L 

=  reach  length. 

V 

=  velocity. 

SLP0 

=  normal  slope. 

D 

=  depth. 

Since  T2  and  C2  are  dependent  upon  02,  an 
iterative  technique  is  required  to  solve  the 
routing  equations.  In  equation  7,  Ia  and  Oi  are 
known,  and  C\  can  be  computed  from  equation 
9.  This  leaves  only  02  and  C2  as  unknowns.  Ox 
can  be  used  as  a  first  approximation  of  O  2 .  The 
normal  depth  and  velocity  for  the  approximate 
value  of  02  are  entered  into  equation  11  for 
computing  T2.  Then  equation  8  is  used  to 
compute  C2.  The  second  approximation  of  02  is 
then  obtained  from  equation  7.  This  iterative 
process  continues  until  the  difference  between 
successive  02  values  is  acceptable.  HYMO  is  set 
to  accept  differences  of  0.1  percent  or  less. 
Usually  about  four  iterations  are  required. 
Reservoirs 

HYMO  uses  the  storage-indication  method  (3) 
to  route  floods  through  reservoirs.  This  method 
has  been  widely  used  and  accepted  because  it  is 
practical  and  accurate.  The  SCS  National  Engi- 
neering Handbook  gives  detailed  instructions  for 
using  the  method. 


Rating  Curves 

Rating  curves  must  be  available  at  enough 
locations  along  a  valley  to  adequately  describe 
the  hydraulics  of  the  stream  and  valley.  Most  of 
these  rating  curves  must  be  computed  because 
there  are  never  enough  measured  rating  curves. 

HYMO  uses  Manning's  equation  to  compute 
the  normal  flow-rating  curves  that  are  used  in 
the  VSC  flood-routing  method.  The  normal 
flood-plain  slope  is  determined  for  each  valley 
section  by  plotting  a  profile  of  the  flood  plain. 
The  normal  channel  slope  is  determined  by 
plotting  a  profile  of  the  flood  plain  with  channel 
distances. 

Sediment  Yields 

The  universal  soil  loss  equation  (8)  was 
modified  to  compute  the  sediment  yield  for 
individual  storms  on  watersheds.  The  modified 
equation  is 

S=95(qp  X<2)0-56  xKxCxPxLS,  (12) 

where     Si  =  sediment  yield  in  tons, 

#n  =  peak  flow  rate  in  cubic  feet  per  second, 
Q  =  volume  of  runoff  in  acre-feet, 
K\  =  the  soil-erodibility  factor, 
C  =  the  cropping-management  factor, 
P  =  the  erosion  control  practice  factor, 

and      LS  =  the  slope  length  and  gradient  factor. 

Detailed  instructions  for  determining  K,  C,  P, 
and  LS  are  given  by  Wischmeier  and  Smith  (8). 

Since  equation  12  was  designed  to  compute 
sediment  yield  from  watersheds,  a  delivery  ratio 
is  not  needed.  The  delivery  ratio  is  built  into 
equation  12  by  including  the  peak  flow  rate. 
Many  of  the  watershed  characteristics  that  in- 
fluence the  peak  flow  rate  also  affect  the 
delivery  ratio.  Equation  12  has  performed  well 
under  limited  testing,  but  future  refinements  are 
expected. 


RULES  FOR  USING  HYMO 


The  reader  should  refer  to  the  table,  "Ex- 
ample Input  for  HYMO  Commands,"  as  he 
follows  the  narrative  description  of  the  rules. 
The  example  HYMO  program  that  is  presented 
near  the  end  of  the  manual  may  also  be  helpful. 

General  Rules 

HYMO  commands  are  expressed  in  the  first 
20  columns  of  the  data  card,  and  columns  21 


through  79  are  used  for  numeric  data  and 
keywords.  Column  80  is  reserved  for  a  page- 
change  code  (an  asterisk  in  column  80  causes  the 
card  to  be  printed  on  a  new  page).  Continuation 
cards  are  allowed  when  59  characters  are  insuf- 
ficient to  express  the  data. 

The  data  can  be  written  in  any  format,  but  at 
least  one  blank  space  must  be  left  between  data 
items.  A  decimal  is  required  for  numbers  con- 
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taining  fractions,  but  not  for  whole  numbers. 
Keywords  can  be  written  with  the  data  to 
describe  individual  data  items.  Comment  cards 


may  be  used  at  any  point  in  a  HYMO  program 
by  punching  an  asterisk  in  column  1  and  the 
comment  in  columns  2  through  79. 


Command 


START 
STORE  HYD 

COMPUTE  HYD 


Example  input  for  HYMO  commands 


Required  input 


RAINFALL  BEGINS  AT  12.5  HRS         PUNCH  CODE=l 

ID=1    HYDNO=301    DT=.2  HR   DA=1.5  SQ  MI 

FLOW  RATES=  0  10  50  100  500  1000  1800  2000  1900  1500  1200 

1000  800  600  500  400  300  200  100  50  10  1 

ID=2   HYD  NO=302   DT=. 5  HR   DA=2.1  SQ  MI  CN=90 
HT=100  FT  L=3.3  MI 

MASS  RAINFALL  =  0  .31  .61  1.04  1.84  2.74  3.06  3.45  4.33  4.75 


PRINT  HYD 

ID=2 

CODE=l 

PUNCH  HYD 

ID=2 

PLOT  HYD 

ID  1=3 

ID  11=4 

ADD  HYD 

ID=4 

HYD  NO=101    ID  1=5   ID  11=6 

STORE  RATING  CURVE 

ID=2 

VS  NO=15 

ELEV 

AREA 

FLOW 

496.6 

0 

0 

497 

2 

1 

498 

9 

19 

499 

19 

52 

500 

30 

98 

COMPUTE  RATING  CURVE 

ID=1  VS  NO=10  NO  SEGS=3  MIN  ELEV=482  FT  MAX  ELEV=492  FT 

CH  SLP=.006   FP  SLP= 

=  .0075   N=.05  DIST=175  FT 

N=-.03   DIST=205  FT 

N=.05  DIST=450  FT 

DIST 

ELEV  DIST 

ELEV     DIST  ELEV 

DIST 

ELEV 

0 

492.0  100 

490.0       175  484.0 

188 

482.0 

190 

482.0  205 

484.0       250  486.0 

275 

488.0 

310 

490.0  450 

492.0 

COMPUTE  TRAVEL  TIME 

ID=3 

REACH  NO=8 

NO  VS=5   L=4500  FT  SLP= 

=.0075 

ROUTE 

ID=3 

HYD  NO=8   INFLOW  ID=6   DT=.25  HR 

ROUTE  RESERVOIR 

ID=5 

HYD  NO=501 

INFLOW  ID=1 

OUTFLOW  (CFS) 

STORAGE  (AC  FT) 

0 

0 

22 

533 

200 

555 

1000 

601 

2000 

648 

3000 

694 

ERROR  ANALYSIS 

ID  1= 

3    ID  11=5 

SEDIMENT  YIELD 
FINISH 


ID=5  SOIL=.34  CROP=.5  EP=.6  LS=.3 
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Six  hydrographs  can  be  stored  in  a  HYMO 
program  at  a  time.  The  hydrographs  are  identi- 
fied by  storage  location  numbers  1  through  6. 
Therefore,  the  same  storage  location  number 
must  be  used  for  many  hydrographs  in  a  HYMO 
program.  This  is  especially  true  when  routing  is 
done  through  large  watersheds.  However,  no 
more  than  six  hydrographs  are  ever  needed  at 
one  time  because  HYMO  programs  begin  at  the 
head  of  a  watershed  and  work  downstream 
through  one  reach  at  a  time.  When  a  storage 
location  number  is  used  to  store  or  compute 
another  hydrograph,  the  first  hydrograph  is  lost. 
The  user  should  be  sure  that  the  hydrograph  will 
not  be  referred  to  again  before  using  the  storage 
location  number  for  another  command. 

To  store,  compute,  or  route  a  hydrograph,  the 
user  must  specify  the  time  increment.  There  are 
no  rigid  rules  about  selecting  the  time  incre- 
ment, but  generally  it  should  not  be  greater  than 
one-fifth  of  the  time  to  the  peak  of  the 
hydrograph.  This  rule  usually  provides  enough 
points  to  adequately  define  the  hydrograph.  All 
hydrographs  are  limited  to  300  points. 

For  the  commands  "STORE  HYD,"  "COM- 
PUTE HYD,"  "ADD  HYD,"  "ROUTE,"  and 
"ROUTE  RESERVOIR,"  the  user  must  specify 
the  number  of  the  outflow  hydrograph.  The 
hydrograph  identification  numbers  are  used  to 
designate  specific  routing  reaches,  incremental 
areas,  reservoirs,  and  partial  hydrographs.  The 
partial  hydrograph  number  is  given  to  all  hydro- 
graphs  other  than  outflow  hydrographs  from 
reaches,  incremental  areas,  or  reservoirs.  The 
identification  numbers  for  each  group  are 

Reaches   1-100 

Partial  hydrographs   101-300 

Incremental  areas   301-500 

Reservoirs   501+ 


Command  Rules 

The  first  command  for  any  watershed  is 
START.  The  two  data  items  associated  with  this 
command  are  the  time  rainfall  begins  on  the 
watershed  and  a  code  for  punching  output  data. 
If  a  storm  is  to  be  routed  through  a  watershed 
only  once,  the  punch  code  is  deleted.  However, 
if  more  than  one  routing  is  to  be  performed,  set 
the  punch  code  equal  to  a  positive  number,  and 
the  output  data  for  the  first  routing  will  be 


punched  for  use  in  the  second  routing.  More 
than  one  routing  is  usually  required. 

Two  commands,  RECALL  HYD  and  STORE 
TRAVEL  TIME,  were  designed  to  be  computer 
punched  for  second  routings;  consequently, 
these  commands  do  not  appear  in  the  table. 

The  STORE  HYD  command  is  used  to  store 
the  coordinates  of  a  hydrograph  in  the  com- 
puter. It  can  be  used  for  storing  measured 
hydrographs  or  hydrographs  computed  by 
methods  other  than  the  one  used  in  HYMO.  The 
input  data  required  for  STORE  HYD  are  storage 
location  number,  hydrograph  identification 
number,  time  increment,  watershed  area,  and 
flow  rates  of  the  hydrograph  at  the  specified 
time  increment. 

The  COMPUTE  HYD  command  is  used  to 
compute  hydrographs  from  the  incremental 
areas  of  the  watershed.  The  first  five  items  of 
data  are  storage  location  number,  hydrograph 
identification  number,  time  increment,  water- 
shed area,  and  SCS  runoff  curve  number  (3). 
Normally,  data  items  6  and  7  are  watershed 
height  and  main  stem  length.  The  height  and 
length  are  used  to  compute  the  recession  con- 
stant K  and  the  time  to  peak  tp  .  However,  if  K 
and  tp  are  known  or  estimated  by  some  other 
method,  they  can  be  entered  directly  into  the 
program.  This  is  accomplished  by  placing  a 
minus  sign  before  the  values  of  K  and  tp  and 
entering  them  as  data  items  6  and  7,  respective- 
ly. The  remaining  data  items  are  values  of  the 
mass  rainfall  at  the  specified  time  increment. 

Since  most  watersheds  have  a  limited  number 
of  rain  gages,  the  same  mass  rainfall  data  may  be 
used  to  develop  several  hydrographs.  Once  the 
mass  rainfall  data  have  been  entered  in  a 
COMPUTE  HYD  command,  they  can  be  re- 
peated for  any  number  of  COMPUTE  HYD 
commands  without  re  punching  the  data. 
Instead,  punch  a  negative  number  for  the  eighth 
data  item  of  all  COMPUTE  HYD  commands  that 
use  the  same  rain  gage.  When  data  from  another 
rain  gage  are  entered,  the  data  from  the  first  rain 
gage  are  lost  and  cannot  be  recalled  by  using  the 
negative  number  code. 

The  RECALL  HYD  command  is  one  of  the 
two  commands  that  are  computer  punched. 
When  the  punch  code  is  a  positive  number,  the 
output  from  STORE  HYD  and  COMPUTE  HYD 
are  punched  on  cards  with  the  RECALL  HYD 
command.  The  RECALL  HYD  command  stores 
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the  computed  and  stored  hydrographs  on  cards; 
it  is  therefore  not  necessary  to  recompute 
hydrographs  for  future  routings.  Instead,  the 
previously  computed  hydrographs  are  read  into 
the  program,  thus  saving  considerable  computer 
time. 

Although  the  input  data  for  the  RECALL 
HYD  command  are  never  punched  manually,  a 
list  of  the  data  items  may  be  helpful  in  checking 
computer-punched  cards.  The  input  data  are 
storage  location  number,  hydrograph  identifica- 
tion number,  time  increment,  drainage  area, 
peak  flow  rate,  runoff  volume,  number  of 
hydrograph  points,  and  flow  rates  of  the  hydro- 
graph. 

The  PRINT  HYD  command  is  used  to  print 
coordinates  of  a  hydrograph,  volume  of  runoff, 
and  peak  flow  rate .  The  required  input  data  are 
the  storage  location  number  and  a  peak-volume 
code.  The  peak-volume  code  is  deleted  if  a 
complete  hydrograph  printout  is  desired.  If  a 
printout  of  only  the  runoff  volume  and  the  peak 
flow  rate  is  needed,  the  peak-volume  code  is  set 
to  a  positive  value. 

The  PUNCH  HYD  command  is  used  to  punch 
iany  hydrograph  in  a  HYMO  program  in  the 
proper  form  for  the  RECALL  HYD  command. 
PUNCH  HYD  has  two  purposes:  (1)  If  the 
punch  code  is  not  used,  PUNCH  HYD  can  be 
used  to  punch  one  or  more  hydrographs  for 
future  use;  and  (2)  if  it  is  desirable  to  punch 
outflow  hydrographs  associated  with  ROUTE, 
ROUTE  RESERVOIR,  or  ADD  HYD,  PUNCH 
HYD  must  be  used  because  the  punch  code  only 
provides  for  punching  hydrographs  associated 
with  STORE  HYD  and  COMPUTE  HYD.  The 
only  datum  required  for  PUNCH  HYD  is  the 
storage  location  number  of  the  hydrograph  to  be 
punched. 

The  PLOT  HYD  command  is  used  to  plot 
hydrographs  in  a  HYMO  program.  It  will  plot 
one  hydrograph  on  a  set  of  axes,  or  if  a 
comparison  is  desired,  it  will  plot  two  hydro- 
graphs  on  the  same  set  of  axes.  The  required 
input  data  are  the  storage  location  numbers  of 
the  hydrographs  to  be  plotted. 

The  ADD  HYD  command  adds  the  coordi- 
nates of  any  two  hydrographs.  The  hydrographs 
are  added  at  a  time  increment  equal  to  that  of 
the  hydrograph  with  the  shorter  time  increment. 
The  only  data  required  are  the  storage  location 
number  and  hydrograph  identification  number 


of  the  added  hydrograph  and  the  storage  loca- 
tion numbers  of  the  two  hydrographs  to  be 
added. 

The  STORE  RATING  CURVE  command  is 
used  to  store  rating  curves  that  have  been 
measured  or  computed  previously.  STORE 
RATING  CURVE  will  save  considerable  com- 
puter time  if  measured  or  computed  rating 
curves  are  available.  The  input  data  are  the 
storage  location  number,  valley,  section  number, 
and  individual  rating  curve  points  described  by 
elevation,  end-area,  and  flow  rate.  The  number 
of  points  used  to  describe  a  rating  curve  is 
limited  to  20. 

The  COMPUTE  RATING  CURVE  command 
is  used  to  compute  the  stage-area-flow  relation- 
ship for  a  valley  section.  The  input  data  are 
storage  location  number,  valley  section  number, 
number  of  segments  in  the  vaHey  section, 
minimum  elevation,  maximum  elevation,  chan- 
nel and  flood -plain  slopes,  Manning's  n  value  and 
segment  boundary  point  for  each  segment,  and 
horizontal  and  vertical  position  of  points  de- 
scribing the  valley  section. 

The  storage  location  numbers  of  the  valley 
sections  in  a  particular  reach  must  begin  with  1 
and  increase  by  one  for  each  valley  section  in 
the  reach.  However,  the  numbers  are  assigned 
without  regard  to  upstream  or  downstream 
order.  The  valley  section  identification  number 
can  be  any  number  from  0.1  to  999.9.  These 
rules  concerning  storage  location  and  valley 
section  identification  numbers  also  apply  to  the 
STORE  RATING  CURVE  command. 

Normally,  valley  sections  are  divided  into 
three  segments  (two  flood-plain  segments  and  a 
channel  segment)  for  computing  the  rating 
curve.  However,  some  valley  sections  may  have 
more  than  one  channel  or  may  have  an  extreme 
variation  in  n  values  across  the  flood  plain,  thus 
requiring  more  than  three  segments.  A  maxi- 
mum of  six  segments  is  permitted.  Manning's  n 
values  for  each  segment  are  input  with  seg- 
ment boundary  point  (distance  from  the  begin- 
ning of  the  valley  section  to  the  end  of  the 
segment).  Flood-plain  n  values  are  positive  and 
channel  n  values  are  negative. 

Twenty  points  are  used  to  define  a  rating 
curve.  The  location  of  the  points  is  determined 
by  dividing  the  difference  between  the  maxi- 
mum and  minimum  elevations  into  19  equal 
increments. 
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The  COMPUTE  TRAVEL  TIME  command  is 
used  to  compute  the  normal  flow  travel  time 
relationship  used  in  ROUTE .  The  input  data  are 
storage  location  number,  reach  identification 
number,  number  of  valley  sections  in  the  reach, 
reach  length,  and  slope.  The  reach  identification 
number  can  be  any  number  from  0.1  to  999.9. 
The  maximum  number  of  valley  sections  per 
reach  is  six.  The  slope  can  be  either  the  channel 
or  flood-plain  slope  or  a  weighted  average  of  the 
two.  If  flow  is  confined  to  the  channel,  the 
channel  slope  is  of  course  applicable.  If  most  of 
the  flow  is  in  the  flood  plain,  usually  the 
flood-plain  slope  is  used.  However,  a  weighted 
slope  based  on  the  relative  rates  of  flow  in  the 
channel  and  the  flood  plain  may  be  used. 

The  COMPUTE  TRAVEL  TIME  command 
considers  each  rating  curve  in  the  reach  in 
computing  the  travel  time  flow  relationship. 
COMPUTE  TRAVEL  TIME  automatically 
selects  the  flow  rates  that  are  used  in  computing 
individual  travel  times.  The  flow  rates  of  the 
rating  curve  with  the  lowest  maximum  flow  rate 
are  chosen.  If  the  flow  rates  of  any  other  rating 
curve  in  the  reach  were  chosen,  the  rating  curve 
with  the  lowest,  maximum  flow  rate  would  have 
to  be  extrapolated.  The  travel  time  table  is 
limited  to  19  points  because  of  the  20-point 
limit  for  rating  curves. 

The  STORE  TRAVEL  TIME  command  is  one 
of  the  two  computer-punched  commands.  When 
the  punch  code  is  a  positive  number,  the  output 
from  COMPUTE  TRAVEL  TIME  is  punched  on 
cards  with  the  STORE  TRAVEL  TIME  com- 
mand. Therefore,  it  is  not  necessary  to  re- 
compute rating  curves  or  travel  time  for  future 
routings.  Instead,  STORE  TRAVEL  TIME  reads 
the  previously  computed  travel  time  flow  rela- 
tionship into  the  program,  thus  saving  consider- 
able computer  time. 

The  input  data  for  STORE  TRAVEL  TIME 
are  not  punched  manually,  but  a  list  of  data 
items  may  be  helpful  in  checking  computer- 
punched  cards.  The  input  data  are  storage 
location  number,  reach  identification  number, 


reach  length,  slope,  and  individual  points  of  the 
relationship  defined  by  depth,  flow,  and  travel 
time. 

The  ROUTE  command  is  used  to  route  floods 
through  streams  and  valleys.  The  input  data  are 
storage  location  number  and  hydrograph  identi- 
fication number  of  the  outflow  hydrograph, 
storage  location  number  of  the  inflow  hydro- 
graph,  and  time  increment.  The  storage  location 
number  of  the  outflow  hydrograph  must  be  the 
same  as  the  storage  location  number  used  in 
COMPUTE  TRAVEL  TIME  for  the  reach.  To 
prevent  unnecessary  program  stoppage,  ROUTE 
extrapolates  the  travel-time  table  when  it  is 
exceeded  and  writes  the  message,  "TRAVEL 
TIME  TABLE  EXCEEDED." 

The  ROUTE  RESERVOIR  command  is  used 
to  route  floods  through  reservoirs.  The  input 
data  are  storage  location  number  and  hydro- 
graph  identification  number  of  the  outflow 
hydrograph,  storage  location  number  of  the 
inflow  hydrograph,  and  individual  points  of  the 
reservoir's  outflow-storage  relationship.  The  out- 
flow-storage relationship  must  be  expressed  in 
20  points  or  less.  If  the  outflow-storage  relation- 
ship is  exceeded,  ROUTE  RESERVOIR  will 
extrapolate  the  relationship  and  write  the  mes- 
sage, "STORAGE-DISCHARGE  TABLE  EX- 
CEEDED." 

The  ERROR  ANALYSIS  command  is  used  to 
determine  the  error  standard  deviation  and  the 
percentage  error  in  peak  flow  between  any  two 
hydrographs  in  a  HYMO  program.  These  func- 
tions make  ERROR  ANALYSIS  useful  in  re- 
search. The  input  data  are  the  storage  location 
numbers  of  the  two  hydrographs  to  be  analyzed. 

The  SEDIMENT  YIELD  command  is  used  to 
compute  the  sediment  yield  at  any  point  in  a 
watershed.  Input  data  required  are  storage  loca- 
tion number  of  the  hydrograph  from  the  area,  a 
soils  factor,  a  crop  factor,  a  slope  length  and 
gradient  factor,  and  a  conservation  practice 
factor  (8). 

The  FINISH  command  is  used  to  end  HYMO 
programs.  There  are  no  data  associated  with 
FINISH. 


EXAMPLE  HYMO  PROGRAM 

A  short  example  problem  is  presented  to  Riesel,  Tex.  A  flood  will  be  routed  through  the 
demonstrate  HYMO.  Figure  4  is  a  map  of  the  watershed  in  its  present  condition,  and  the 
6.84-square-mile  Brushy  Creek  watershed  near    routed  outflow  hydrograph  will  be  compared  to 


3 


the  hydrograph  measured  at  gaging  station  G. 
Also  the  sediment  yield  will  be  predicted  and 
compared  with  the  measured  sediment  yield. 
Then  the  same  flood  will  be  routed  through  the 
watershed  with  two  proposed  reservoirs.  To 
determine  the  effects  of  the  reservoirs,  the 
outflow  hydrograph  and  sediment  yield  will  be 
compared  to  the  outflow  hydrograph  and  sedi- 
ment yield  of  the  present-condition  routing. 


Comment  cards  and  keywords  are  used  liber- 
ally in  the  example  problem  to  acquaint  the  user 
with  HYMO.  After  becoming  familiar  with 
HYMO,  the  user  may  write  fewer  comments  and 
keywords,  but  generally  users  find  them  both 
quite  helpful  in  describing  the  problem.  To  save 
space  in  the  example  problem,  few  of  the 
hydrographs  are  printed  or  plotted.  Some  users 
may  choose  to  print  and  plot  all  hydrographs. 


SCALE 
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Figure  4.  Brushy  Creek  watershed  near  Riesel,  Tex. 
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APPENDIX 

HYMO  MAIN  PROGRAM  AND  SUBROUTINES 


//SOPTIONS 


C 

HYMO 

A 

1 

c 

A 

2 

c 

THIS  PROGRAM   IS  A  PROBLEM  ORIENTED  COMPUTER  LANGUAGE  FOR 

A 

3 

c 

DEVELOPING  WATERSHED  MODELS. 

A 

4 

c 

A 

5 

c 

A 

6 

c 

DEFINITION  OF  VARIABLES 

A 

7 

c 

ID  =  STORAGE  LOCATION  NUMBER 

A 

8 

c 

NHD  =  HYDROGRAPH   IDENTIFICATION  NUMBER 

A 

9 

\ 

c 

DT   =  TIME    INCREMENT    IN  HOURS 

A 

10 

f 

c 

DA    =   DRAINAGE    AREA    IN    SQ  MI 

A 

1 1 

c 

OCFS  =  FLOW  RATE  ARRAY 

A 

12 

c 

□  CFSF  =  CUTOFF  FLOW   RATE    IN  CFS 

A 

13 

c 

trim    _     k  1 1  i  m  n  r  n     fi  t     n  n  •*  kt  T  c     Vki     a     i  i  u  r\n  nr  n>  a  niii 

I  END  =   NUMBER  OF  POINTS   IN  A  HYDROGRAPH 

A 

14 

c 

C  f  C~  C~               r(   nu     n  A  T  r      A  D  O  A  \S     C  n  □      TD  A  l(  C  1       TTUC     T  A  □  t  C 

SC  rS  =    FLOW  R Alt    ARRAY   rOK    TRAVEL   TIME  TABLE 

A 

15 

c 

A      _     run     ft  r%  c  A      i  \i     c/\  CT 

A    =  END  AREA    IN   5Q  FT 

A 

16 

c 

Q   =  FLOW  RATE   ARRAY  FOR  RATING  CURVES 

A 

17 

c 

RAIN  =  MASS  RAINFALL  ARRAY 

A 

1 8 

1 

CO MMON  CFS(300),OCFS(300,6) ,1  END (6ltDATA(310)fDA(6l,DP(20l,NPU»NHD 

A 

19 

If SCFS(20lfCI20) ,A(20,6),Q(20,6)«RAIN(200) »DEEP(20»6)»NER»MAXN0,NC0 

A 

20 

2MMfICCfNC0DEfDIST(6)fSEGN(6)fCTBLE(50flllfITBLEi50f2lfZALFA(20)fDT 

A 

21 

3( 6)»TIME»PEAK(6)»R0IN»ISG(6) 

A 

22 

2 

NC  CDE=0 

A 

23 

c 

READ  ZALFA  ARRAY 

A 

24 

3 

READ  I5f20)    1  ZALF  A(  I )  f  1=  I  f  1  5) 

A 

25 

c 

READ  NUMBER  OF  COMMANDS 

A 

26 

4 

READ   (5,21 )  NCOMM 

A 

27 

c 

READ  COMMAND  TABLE 

A 

28 

5 

n>—  Ar%      1   r       ~\  ~\  ft        */^Trtir*#t         1  ft         ■       1        «  1    ft       *    t  t  r»  *  J  v         ■  ft         ■       «       ~%  \       *         1      ais*  <~ft  aa  >  a  ft 

READ  (  5,22)    (  (CTBLE  (  I  ,J),J=1,  11)  ,(  ITBITtl  I,J),J=1,2),I  =  1,NC0MMI 

A 

29 

6 

1   1  fi     IT  "S-  j—         fly          "*    H    6          4    ~*    A  i     «~  A    J   "aT    1         T         V          11    C  4 

WRITE   (6,23)    I ZALF A ( I ) , I =1 , 15) 

A 

30 

7 

WR  ITE  (6,24) 

A 

31 

8 

i  1  r,    17  ir  r-        #    v          1  r  t         f    1  /*  T  n  1     f*"i    T          aft          a        S         ti%         it  T  ri  ■    r-  1  »           1  w          ■       «        Oft        77        1        ii«>       aa  aa  ft 

WRITE    (6,25)    ( ( CTBL E 1  I , J )  , J  =  l , 1 1 ) , I ITBLE ( I  , J ) , J  =  lf 2 ) , I  =  1 , NCOMM ) 

A 

32 

c 

ZERO  CODES   —   NPU  =   PUNCHf    ICC   =   CONTINUATION   CARDf    NER  =  ERROR 

A 

33 

9 

NPU=0 

A 

34 

1 0 

I  C  C=0 

A 

35 

1 1 

1 

N  E  R=0 

A 

36 

12 

CALL  HONDO 

A 

37 

1 3 

t  j—       1  ur  n  ft       ^       0  ft 

IF    (NER)  2f2fl9 

A 

38 

14 

2 

GO   TO   1 3  f  4  ,  5  , 6  ,  7,  8  ,9  ,  10  ,  1 1 ,  12,  13  , 1  4, 15  , 16,  17,  1 8  ,  19  )  ,  NCODE 

A 

39 

1  5 

3 

Tffftftf        fiATAfft  ft 

T I ME=DATA ( 1 ) 

A 

40 

16 

NPU=DATA( 2 ) 

A 

41 

17 

GO  TO  1 

A 

42 

1  8 

4 

CALL  STHYD 

A 

43 

19 

GO  TO  1 

A 

44 

20 

5 

CALL  RECHD 

A 

45 

2 1 

GO   TO  1 

A 

46 

22 

6 

CALL  CMPHYD 

A 

A 

4  1 

23 

GO   TO  1 

A 

A 

y.  q 
•fO 

24 

7 

CALL  PRTHYD 

A 

A 

/.  y* 

49 

25 

GO  TO  1 

A 

A 

30 

26 

8 

CALL  PUHYD 

A 

A 

5 1 

27 

GO  TO  1 

A 

A 

c  ■> 

D<: 

28 

9 

CALL  HPLOT 

A 

53 

29 

GO  TO  1 

A 

54 

30 

10 

CALL  ADHYD 

A 

55 

31 

GO   TO  1 

A 

56 

32 

11 

CALL  SRC 

A 

57 

33 

GO   TO  1 

A 

58 

34 

12 

CALL  CMPRC 

A 

59 

54 


-n  ^-i  c  n  C  73 
§•  £  </>  ^  2.  5, 


\ 


35 

GO  TO  1 

A 

60 

36 

13 

CALL  STT 

A 

6 1 

37 

GO  TO  1 

A 

62 

38 

14 

("•11            f    Ik  f\  T  T 

CALL  >CMPTT 

A 

63 

39 

GO  TO  1 

A 

64 

40 

1  5 

CALL  ROUTE 

A 

65 

41 

GO   TO  1 

A 

66 

42 

16 

CALL  RESVO 

A 

6  7 

43 

GO  TO  1 

A 

6  8 

44 

17 

CALL  ERROR 

A 

69 

45 

GO  TO  1 

A 

70 

46 

18 

CALL  SEDT 

A 

71 

47 

GO   TO  1 

A 

72 

48 

19 

STOP 

A 

73 

C 

A 

74 

49 

20 

FORMAT  (15A1) 

A 

75 

50 

2  1 

FORMAT  (12) 

A 

76 

5  1 

22 

FORMAT    ( 2 Al , 9A2 ,2 13 ) 

A 

77 

52 

23 

FORMAT    ( 1H1 ,9X, 8HZALFA  =  ,15A1///) 

A 

78 

53 

24 

FORMAT    { 16X,13HC0MMAND  TABLE//) 

A 

79 

54 

25 

FORMAT    ( 10X  ,2A 1 ,9A 2 ,2  I  3 ) 

A 

80 

55 

END 

A 

8  l- 

56 

SUBROUTINE  HONDO 

B 

1 

C 

THIS  SUBROUTINE  READS    IN  A  DATA  CARD   ,    SEARCHES  AN  ALPHAMERIC 

B 

2 

C 

CODE  TABLE   TO  DETERMINE   THE   NCODE  OF  THE  OPERATION  AND 

B 

3 

c 

COLLECTS  VARIABLES  FROM  THE  FREEFLOAT ING  DATA  FIELD 

B 

4 

5  7 

COMMON  CF  S ( 300 ) ,OCFS( 300,6) ,IEND(6),DATA(310),DA(6),DP(20),NPU,NHD 

B 

5 

If  SCFS ( 201 1 C ( 20 ) t A( 20, 61 , Q ( 20, 61 , RAIN (200) , DEEP (20»  6) tNER , MAXNO, NCO 

B 

6 

2MM,ICC,NC00E,DIST(6),SEGN(6),CTBLE(50»11),ITBLE(50,2),ZALFA(20)»DT 

B 

7 

3(6),TIME,PEAK(6),R0IN,ISG(6) 

3 

8 

58 

DIMENSION  CHAR ( 60 ) ,  ALPHA(ll) 

B 

9 

59 

DIMENSION  AUXA(IO),  AUXB(IO) 

B 

10 

60 

I F    (  ICC )  1,1,3 

B 

1 1 

c 

READ  IN  DATA  CARD 

B 

12 

6  1 

1 

READ  (5,42)    (  ALPHA( I  I ,1  =  1,11) ,(CHAR( I ) ,1  =  1, 60) 

B 

1  3 

c 

IF   FIRST  CHARACTER    IS   BLANK  THE  CARD   IS  A  CONTINUATION  OF 

B 

14 

c 

PREVIOUS  CARD. 

B 

1  5 

62 

IF    ( AL PHA(  1  )-ZALF A(  11 ) )  2,9,2 

B 

16 

63 

2 

IF    (  ICC  )   3 , 3,40 

B 

17 

C 

ASTERISK   IN  COL.    80   MEANS   SKIP  TO  NEW  PAGE   BEFORE   PRINTING  CARD 

B 

18 

64 

3 

IF    ( CHAR ( 60 )-ZALF A( I 1 ) )  4,5,4 

B 

19 

65 

4 

WRITE   ( 6,43 1 

B 

20 

66 

5 

WRITE    (6,44)    ( ALPHA ( I ) , 1= 1 , 1 1 ) , ( CHAR  (  I  )  ,  1  =  I  ,60 > 

B 

21 

C 

IF  FIRST  CHARACTER    ISA*     THE   PREVIOUS  CARD  WAS  A   COMMENT  CARD 

B 

22 

67 

IF    ( ALPHA( 1 ) -ZALF A( 12 ) )  10,6,10 

B 

23 

C 

IF    PUNCH  CODE   POSITIVE,   COMMENT   CARDS   ARE  PUNCHED. 

B 

24 

68 

6 

IF    (NPU)  8,8,7 

B 

25 

69 

7 

WRITE    (7,45)    ( ALPHA ( I ) , 1= I , 11 ) , ( CHAR ( I ) , 1=1 ,60 ) 

B 

26 

70 

8 

ICC=0 

B 

27 

71 

GO  TO  1 

B 

28 

72 

9 

WRITE    (6,44)    ( ALPHA( I ) , 1=1 , 1 1 ) , ( CHAR ( I ) , 1= 1 ,60 ) 

B 

29 

73 

GO  TO  24 

B 

30 

C 

SEARCH  FIRST  TWO  ALPHAMERIC   CHARACTERS  TO   SEE    IF  THEY  ARE  NUMBERS 

B 

31 

74 

10 

ICC=1 

B 

32 

75 

DO   12  1=1,10 

B 

33 

76 

IF    ( ALPHA(  l)-ZALFA(  I)  )  11,15,11 

B 

34 

77 

1 1 

IF    ( ALPHA( 2)-ZALFA( I ) )  12,15,12 

B 

35 

78 

12 

CONTINUE 

B 

36 

C 

STATEMENT    NUMBER    7    IS    BRANCHED   TO    IF   NUMBERS   ARE  PRESENT 

B 

37 

55 


c 

IF   NOT   NUMBER  SEARCH  COMMAND  TABLE 

FOR  MATCH 

B 

J  o 

c 

CALL  FIRST    10  VALUES  FROM  PERMANENT 

DATA 

STORAGE 

g 

39 

79 

DO   14  1=1,50 

g 

80 

DO   13  J  =  l,  11 

g 

41 

8  1 

c 

IF   (CTBLE(  I  , J  )-ALPHA( J) )  14,13, 

SN   10=PART  MATCH 

14 

B 

42 
43 

82 

1  3 

CONTINUE 

g 

44 

c 

IF   THIS  LOOP    IS  COMPLETED  WE  HAVE 

COMPLETE 

MATCH- 

CALL  NCODE 

B 

45 

c 

AND  MAX  NUMBER   AND  EXIT  LOOP 

g 

46 

83 

NC  ODE= I TBL  E ( I ,1) 

g 

47 

84 

MAXNO= ITBLEI I ,2) 

B 

48 

85 

GO  TO  21 

g 

49 

86 

14 

CONTINUE 

g 

50 

c 

IF   MAJOR  LOOPS  FINISHED  WITHOUT 

A 

MATCH 

WRITE 

ERROR 

MESSAGE 

g 

51 

c 

ANO  SET  NER     =  1 

B 

52 

87 

NER=1 

g 

53 

88 

WRITE  (6,46) 

B 

54 

89 

RETURN 

B 

55 

c 

CONVERT   DIGIT    INPUT  CODE  FROM  ALPHAMERIC 

TO 

INTEGER 

FORM 

B 

56 

90 

15 

NCODE=GIT<  ALPHA  ,1,2,1. 1*0.5 

B 

57 

C 

FIND  MAX  NUMBER   OF   DATA   ITEMS  FOR 

THIS  NCODE 

B 

58 

91 

DO   17   1=1, 50 

B 

59 

92 

IF   ( ITBLEI  1 ,1 l-NCODE )  17,16,17 

B 

60 

93 

1  6 

MA  XNO= I TBL  E (  I »  2  ) 

B 

6 1 

94 

GO  TO  21 

B 

62 

95 

17 

c 
c 

CONTINUE 

SEARCH  DATA  ROUTINE 

SEE   IF   ANY   DATA  FOR   THIS  CARD 

B 
B 
B 

63 
64 
65 

96 

DO   19  1=1,50 

B 

66 

97 

IF    ( ITBLE( I ,1 J-NCODE)  19,18,19 

B 

67 

98 

1  8 

M A  XNO= I TBL  E ( I  ,2) 

B 

68 

99 

GO  TO  20 

B 

69 

100 

19 

CONTINUE 

B 

70 

101 

20 

CONTINUE 

B 

71 

10  2 

2  1 

IF    (MAXNO)  23,22,23 

B 

72 

10  3 

22 

c 

RE  TURN 

ZERO  ARRAYS   AN C  COUNTERS 

B 
B 

73 
74 

104 

23 

DO  47  1=1,310 

B 

75 

105 

4  7 

DATA  (I  )=0. 

B 

76 

106 

ND ATA= 1 

B 

77 

107 

24 

NC  HAR=0 

B 

78 

108 

25 

DO  26  1=1,10 

B 

79 

109 

AUXA( I )=0. 

B 

80 

110 

26 

AUXB(  I  )  =0  . 

B 

81 

111 

IT  1=1 

B 

82 

112 

I  T  2=1 

B 

83 

113 

SIGN=1. 

B 

84 

114 

LOGIT=0 

B 

85 

115 

KDGI  T=0 

B 

86 

c 

CARRY  OUT  DIGIT   BY  DIGIT  SEARCH 

AND 

ACCUMULATION 

B 

87 

116 

2  7 

NCHAR=NCHAR+1 

B 

88 

HAVE  WE  CONSIDERED  ALL  CHARACTERS 

RETURN 

IF 

SO 

B 

89 

117 

IF    (NCHAR-60)  28,32,1 

B 

90 

118 

28 

DO  29  1=1,15 

B 

91 

119 

IF    (CHAR(NCHAR)-ZALFA(I ))  29,30 

,29 

B 

92 

120 

29 

CONTINUE 

B 

93 

121 

GO  TO  32 

B 

94 

122 

30 

GO   TO  (33,33,33,33,33,33,33,33, 

33, 

33,32, 

27, 

36, 

32, 

31 

.27),  I 

B 

95 

C 

SN  39     HANDLES  SIGN  CONTROL  ON 

1130 

VERSION 

B 

96 

123 

31 

SIGN=-1.0 

B 

97 

56 


124 

GO  TO  27 

B 

98 

C 

CHARACTER    IS   BLANK   OR  COMMA  -   DOES  IT 

FOLLOW 

A  DIGIT 

B 

99 

125 

32 

GO   TO   (27,48),  ITl 

B 

100 

C 

r  ii  ad  at  t  cn     t  c     a    n  i  n  t          ua  f     A     ncr  iuai 
CHAKACicK    IS    A   U  lbl  T    —    HAS    A  UtLlMAL 

BEEN  ENCOUNTERED 

Q 

D 

101 

126 

3  3 

GO   TO    (34,35),  IT2 

B 

102 

127 

34 

LDGI T=LDGI  T  +  l 

B 

103 

128 

I  T  1=2 

B 

104 

129 

AUXA(LOGIT ) =CHAR ( NCHAR ) 

Q 

D 

103 

1  30 

GO   TO  27 

B 

106 

131 

3  5 

KD  G I T=KDG I T* 1 

B 

10  f 

132 

AUXB(KDGIT)=CHAR(  NCHAR  ) 

B 

108 

133 

GO    TO  27 

B 

109 

C 

CHARACTER    IS    A   DECIMAL  -  DOES    IT  FOLLOW  A  DIGIT 

B 

110 

134 

36 

GO   TO   (37,38) ,  ITl 

B 

111 

135 

3  7 

I  Tl  =  2 

B 

112 

136 

L  D  G  I  T=  1 

B 

113 

137 

3  8 

I  T2  =  2 

B 

114 

138 

GO   TO  27 

B 

115 

C 

ROUTINE   TO   CONVERT   ALPHABETIC   ARRAY   TO  FLOATING  POINT 

hi  hi  n  r  Pi 

NUMBER 

B 

116 

139 

4  8 

OA  T A   (NDATA )  =  GIT( AU XA , 1 , LOG  IT, 1 . )  *G I T ( AUX8  , 1, 

10,0. ) 

B 

I  1  f 

140 

DATA  ( NDATA)=DATA(NDATA)*SIGN 

B 

118 

c 

IS  ALL    DATA  FURNISHED     YES-RETURN  NO 

INCREASE  N  DATA 

KEEP  ON 

B 

119 

141 

IF    ( NDA  TA-MAXNO )  41,39,39 

B 

120 

142 

39 

ICC  =  0 

B 

121 

14  3 

40 

RETURN 

B 

122 

144 

4  1 

ND ATA=NDATA*1 

B 

123 

145 

GO   TO  25 

B 

124 

c 

B 

125 

146 

42 

FORMAT    (  2  A 1  ,9  A2  , 60A  1 ) 

B 

126 

147 

4  3 

FORMAT  (1H1) 

B 

127 

148 

44 

FORMAT    (  5X  ,2A1, 9A2, 60A1 ) 

B 

128 

149 

4  0 

FORMAT    (2A1  ,9A2,60A 1) 

B 

129 

1  50 

46 

FORMAT    ( 10X,20HCOMMAND  NOT    IN  TABLE) 

B 

130 

151 

ENC 

B 

131- 

152 

FUNCTION  GIT    ( TC ARD  ,  J , JLA ST , SHI  FT ) 

C 

1 

1 53 

DIMENSION  TCARD(IO),  A(10> 

C 

2 

1 54 

DATA  A(1)/1H1/,A(2)/1H2/,A(3)/1H3/, A(4)/1H4/, 

A(5)/1H5/ 

,  A(6)  /H  6/ 

C 

3 

155 

DATA  A(7)/1H7/,A(8» /l H8 / , A ( 9 ) /1H9/ , A( 10 ) / 1  HO/ 

C 

4 

156 

GI T=0. 

c 

5 

1 57 

TEN=10. 

c 

6 

15  8 

SUM=0. 

c 

7 

1 59 

DO   3  JNOW=J,JLAST 

c 

8 

160 

TTEST=TCARD( JNOW) 

c 

9 

C 

CHECK  FOR  LAST  ENTRY 

c 

10 

161 

IF    { TTEST. EQ.O. )    GO  TO  4 

c 

1 1 

C 

FIND  NUMBER  AND  COMPUTE  VALUE 

c 

12 

162 

DO  2  NUMB=1,10 

c 

1  3 

163 

IF    ( TTEST- A( NUMB ) )  2,1,2 

c 

14 

164 

1 

Z  TEST=NUMB 

c 

1  5 

165 

IF   ( ZTEST.EQ. 10. )  ZTEST=0. 

c 

16 

166 

SUM=SUM*TEN+ZTEST 

c 

17 

GO  TO  3 

c 

1 8 

168 

2 

CONT  INUE 

c 

19 

169 

3 

CONTINUE 

c 

20 

170 

4 

IF    (SHIFT)  6,5,6 

c 

21 

171 

5 

FI =JNOW-l 

c 

22 

172 

SUM=SUM*(0.1**FI  ) 

c 

23 

173 

6 

GI T=SUM 

c 

24 

174 

RETURN 

c 

25 

57 


175  END  C  26- 

176  SUBROUTINE   STHYD  D  ± 
C  THIS  SUBROUTINE   STORES  THE  COORDINATES  OF   HYDROGRAPHS.  D  2 

177  COMMON  CPS (300) ,OCFS( 300,6) , IEND(6) ,DATA(310) ,DA(6> ,DP(20) ,NPU,NHD     D  31 
1, SCFS( 20 ) ,C(20)  ,A(20,6) ,Q(20,6) , RA I N ( 200 ) , DE EP <  20 , 6 ) , NER , MA XNO ,NCO     D  4 
2MM  ,ICC,NC0DE,DIST<6> ,SEGN(6) ,CTBLE (50,11),  I TBL E ( 50 , 2 ) , ZALF A ( 20 ) , DT     D  5 
3(6), TIME, PEAK(6) ,R0IN,ISG(6>  0  6 

178  I D=DAT A ( 1 )  D  7 

179  NHD=DATA( 2 )  D  8 

180  DT(ID)=DATA(3)  0  9 

181  DA( I0)=DATA(4)  D  10 

182  J  =  5  D  u 
C  REMAINING   DATA  ARE   FLOW  RATES  q  12 

183  OCFSd  ,  ID)  =DATA(  J)  0  13 

184  PEAKd  D)  =  l.  D  14 

185  R  0  =DAT  A ( J )  D  15 

186  DO   4   1=2,300  D  16 

187  J=J+1  D  17 

188  OCFS( I , ID) =DATA( J )  D  18 

189  RO=RO+OCFS(I, ID)  D  19 
C  IS  FLOW  RECEDING  0  20 

190  IF  (OCFS( I  ,  IO)-OCFS( 1-1, ID) )  1,2,2  D  21 
C  HAS  FLOW  RECEDED  TO  CUTCFF   RATE  D  22 

191  1  IF  (OCFSd, ID))  5,5,4  D  23 
C  DETERMINE   PEAK  FLOW  D  24 

192  2  IF    (OCFS( I  ,  ID)-PEAK(ID) )   4,4,3  D  25 

193  3  PEAK(ID)=OCFS( I »ID)  D  26 

194  4  CONTINUE  D  27 

195  5  IEND(ID)=I-1  0  28 

196  M  =  I  END (ID)  0  29 

197  RO I N=( RO*DT( ID) ) / ( DA( ID)* 64 5. 333)  D  30 
C  PUNCH  CODE  D  3{ 

198  IF    (NPU)    7,7,6  0  32 

199  6  WRITE   (7,8)    I  0 , NHD , DT ( I D ) , D A ( I D ) , PE AK ( I D ) , ROI N , I END( I D )  D  33 

200  WRITE   (7,9)    ( OCFS( J ,  ID)  ,  J=l , M)  0  34 

201  7  RETURN  D  35 

C  D  36 

202  8  FORMAT!        'RECALL  HYD ' , T2 1 , ' I D= * , 1 1 , T2 9 , ' H YD  N0= * , I  3 » T42 » ' DT=  '  , F9 .     0  37 

16, .»  HRS*  ,T61,  «DA=«  ,  F8.3,  •    SQ  M I  •  /T  2 1 ,  •  PE  AK=  • ,  F7.  0  ,  •  CF  S  «  ,  T  40  , '  R  0=  •  ,     D  38 

2F6.3,'    INCHES • ,T59, "NO  PTS=  '  ,  I  3/ T2 1 , • FLOW   RATES')  D  39 

203  9  FORMAT   (T21,7F8.0)  0  40 

204  ENO  0  41_ 

205  SUBROUTINE  RECHD  E  1 
C  THIS  SUBROUTINE  RECALLS  PREVIOUSLY  COMPUTED  AND  PUNCHED  E  2 
C  HYDROGRAPHS  E  3 

206  COMMON  CFS (300) ,OCFS( 300,6) , IEND(6) ,DATA(310) ,DA(6) ,DP( 20) , NPU, NHD  E  4 
1,SCFS(20),C(20),A(20,6),Q(20,6) , RAIN (200) , DEEP (20, 6) , NER , MA XNO, NCO  E  5 
2MM, ICC,NC0DE,DIST(6) , S EGN ( 6 ) , CT BLE ( 50 , 1 1 ) , I TBL E ( 50 , 2 ) , Z ALF A ( 20 ) , DT  E  6 
3(  6),TIME,PEAK(6) ,ROIN, ISG(6)  E  7 

207  I D=DAT A ( 1 )  E  8 

208  N H  C=DA  T A ( 2 )  E  9 

209  DT ( I D) =DAT A( 3  )  E  10 

210  DA( ID) =DATA(4 )  E  u 

211  PEAKd  D)=DATA(  5)  E  12 

212  R  0 1 N  =  D ATA ( 6 )  E  13 

213  I E  ND( I D)=DATA ( 7)  E  14 

214  M=IEND(IO)  E  15 

215  J=8  E  16 


58 


C 

REMAINING   DATA  ARE   FLOW  RATES 

E 

17 

216 

00    1    1=1, M 

E 

18 

217 

OCFS( I  ,ID)=DATA(J) 

E 

19 

218 

1 

J  =  J+l 

E 

20 

219 

RE  TURN 

E 

21 

220 

END 

E 

2  2- 

221 

SUBROUTINE  CMPHYD 

F 

1 

C 

THIS    PROGRAM   DEVELOPS   A  UNIT  HYDROGRAPH,    CONVERTS   MASS  RAINFALL 

F 

2 

C 

TO   POINT   RUNOFF,    AND  COMPUTES   STORM  HYDROGRAPHS   BY  SUMMATION. 

F 

3 

222 

COMMON  CFS { 300) ,OCFS( 300, 6) ,IEND(6),DATA(310),DA(6),DP(20),NPU,NHD 

F 

1, SCFSI 20) ,C(20)  ,A(20,6) ,Q(20,6) , RAIN (200) , DEEP ( 2C , 6 ) , NER , MA XNO , NCO 

F 

5 

2MM  ,ICC ,NCODE, DI  ST (6) , SEGN (6 ) , CTBLE (50, 1 1 ) , ITBLE( 50,2 ) , Z ALF A ( 20 ) , DT 

F 

6 

3(  6)  ,TIME,  PEAM6)  ,  RO  IN,  ISG(  6) 

F 

7 

223 

I D=DATA ( 1 ) 

F 

8 

224 

NHD=DATA( 2 ) 

F 

9 

225 

D  T ( I D ) =DAT  A ( 3  ) 

F 

10 

226 

DA ( ID)=DATA(4) 

F 

1  1 

227 

CN=DATA( 5  ) 

F 

12 

c 

ARE  K   AND   TP   FURNISHED   OR    WILL    THEY  BE  COMPUTED 

F 

13 

228 

IF    ( DA  T A ( 6 ) )  1,2,2 

F 

14 

229 

1 

XK=-OATA( 6 ) 

F 

15 

230 

TP  =— DAT  A ( 7  ) 

F 

16 

231 

GO  TO  3 

F 

17 

232 

2 

HT=DATA( 6) 

F 

18 

233 

XL  =DATA ( 7  ) 

F 

19 

234 

SLOPE=HT/XL 

F 

20 

235 

XLDW=(  XL**2. ) /DA( ID) 

F 

21 

236 

XK=27.0*(DA(ID)**.231)*(SL0PE**(-.7  77) ) * ( XL  DW**  .124) 

F 

22 

237 

TP=4.6  3*( 0A( ID)**. 4 22 >*<SLOPE**(-. 46) ) * { XL OW** . 1 33 ) 

F 

23 

238 

3 

PEAK! ID) =1. 

F 

24 

239 

DO  4  1=1,300 

F 

25 

240 

4 

OCFS( I , ID) =0. 

F 

26 

C 

COMPUTE   N   BY  ITERATION. 

F 

27 

241 

XN=5.0 

F 

28 

242 

XK  TP=XK/TP 

F 

29 

243 

DO  6  1=1,50 

F 

30 

244 

Tl NF=1 ,*SQRT( 1. /( XN-l . ) ) 

F 

31 

245 

XNl=.05/( XKTP*( ALOG( T I NF/ ( T I NF + .05) )+.05) ) +1. 

F 

32 

246 

DI  FF=ABS(  XM-XN) 

F 

33 

247 

IF    (DIFF-.OOl)  7,7,5 

F 

34 

248 

5 

XN=XN1 

F 

35 

249 

6 

COMTINUF 

F 

36 

250 

WRITE  (6,29) 

F 

37 

251 

GO  TO  28 

F 

38 

C 

DETERMINE  CI. 

F 

39 

252 

7 

DELT=T INF/ 100  . 
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IF    (DCFSM.ID1I  — 02)  28.28.77 
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27 

CFS(K)=ZER0 

I  80 

456 

28 

CONTINUE 

I  81 

457 

WRITE    (6,44)    DOT, (CFS(  I  )  ,1=1, MAX) ,D0T 

I  82 

458 

IF    (ID2)  34,34,29 

I  83 

459 

29 

K  =  l 

I  84 

460 

00   33    1=2, M 
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461 

K=K+  MR  TO 

I 

86 

462 

IF   (OCFS { I , ID2 )-Ql )  30,31,32 

I 

87 

463 

30 

IF   (OCFS( I , ID  2)-Q2)  32,32,31 

I 

88 
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31 

CFS(K)=PLUS 
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GO  TO  33 

I 

90 

466 

32 

CFS( K) =BLANK 

j 

91 

46  7 

33 

CONTINUE 
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46  8 

WRITE    (6,42i    (CFS( I ) , 1=1, MAX) 
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93 

469 

34 

IF    (J-Jl)  36,35,36 

j 
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470 

35 

J 1=J 1+ 10 
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471 

WRITE   (6,43)  02 
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36 

Q1=Q2 

I 

97 

473 

37 

CONTINUE 
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98 

474 

CF  S( 1)=TIME 

j 

99 

475 

DTT=DT(  ID!  )*10./XMRT0 
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100 

C 

PUT  TIME   ARRAY    IN  CFS  AND  WRITE 

TIME  SCALE 

I 

101 

476 

DO  3  8  1=2,12 

I 

102 

477 

38 

CF  S(  I  )=CFS  (  I-D+DTT 

I 

103 

478 

WRITE   (6,45)    (CFS(  I  ),  1  =  1,  12) 

I 

104 

479 

WRITE  (6,46) 

I 

105 

480 

IF    (NPU)  40,40,39 

I 

106 

481 

39 

WRITE   (7,47)  ID1,ID2 

I 

107 

482 

40 

RETURN 

I 

108 

C 

I 

109 

483 

41 

F0RMAT(1X,F6.0,'  CFS.',119A1) 

I 

110 

484 

42 

FORMAT    ( 1H+,  1  IX,  118A1 ) 

I 

111 

485 

43 

FORMAT    (  1H  +  ,F6.0 ) 

112 

486 

44 

FORMAT  (11X,120AII 

[ 

113 

487 

45 

FORMAT  (6X,12F10.3) 

I 

114 

488 

46 

F0RMAT(49X,«TIME  HOURS'///) 

I 

115 

489 

47 

FO  RM  AT (        'PLOT  HYD • ,T21 , '  ID 

1  =  ' 

,11 ,T29,« 

ID  1 1  =  •  ,  1 1 1 

I 

116 

490 

END 
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491 

SUBROUTINE  ADHYD 

J 

1 

C 

THIS  SUBROUTINE  ADDS  TWO  HYDROGRAPHS. 

J 

2 

492 

COMMON  CFS (300)  ,OCFS( 300,6) , I  END (6) , D AT A ( 3 10 ) , D A( 6 ) , DP ( 20 ) , NPU , 

NHD 

J 

3 

1, SCFS( 20) ,C( 20) ,A(20,6> ,Q(20, 

6), 

RAIN(200) 

,DEEP(20,6) , 

NER , MAXNO , 

NCO 

J 

4 

2MM ,1 CC ,NC0DE,DIST(6) , SEGN(6) , 

CTBLE(50,11) 

,  ITBLE(50,2) 

,ZALFA(20) 

,DT 

J 

5 

3( 6),TIME,PEAK(6) ,R0IN,ISG(6) 

J 

6 

493 

I D=DATA( 1) 

J 

7 

494 

NHD=DATA(  2  ) 

J 

8 

495 

I D  1  =  DATA  (  3 ) 

J 

9 

496 

I D2=DATA(4) 

J 

10 

497 

PEAK(ID)=1. 

J 

11 

C 

MAKE  TIME    INCREMENTS   EQUAL  IF 

NOT   EQUAL . 

USE  SMALLER 

INCREMENT 

J 

12 

498 

IF   ( DT( ID1 )-DT( ID2) )  1,3,2 

J 

13 

499 

1 

DT( ID)=DT( IDl ) 

J 

14 

500 

L=  IDl 

J 

15 

50  1 

K=  ID2 

J 

16 

502 

GO  TO  6 

J 

17 

503 

2 

DT ( I D) =DT { ID2 ) 

J 

18 

504 

L=  ID2 

J 

19 

505 

K  =  ID1 

J 

20 

506 

GO  TO  6 

J 

21 

507 

3 

DT( ID)=DT( IDl ) 

J 

22 

508 

IF    ( IEND( ID1)-IEND( ID2) )  4,4, 

5 

J 

23 

509 

4 

M3=IEND( IDl) 

J 

24 

510 

K l=ID2 

J 

25 

511 

IEND(ID)=IEND(ID2> 

J 

26 

512 

GO  TO  18 

J 

27 
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<?13  5  M3  =  IEND(ID2)  J  ^H 

514  K1  =  ID1  J  29 

515  IEND( ID)= I END( ID1)  J  30 

516  GO   TO   18  J  31 
C           DETERMINE   DURATIONS  OF   FLOW  J  32 

517  6  X  I END1= IEND( ID1 )-l  J  33 

518  X  I  END2=IENDU  D2)-l  •  J  34 

519  DUR1=X IEND1*DT( ID1)  J  35 

520  DUR2=XIEND2*DT( I  D2 )  J  36 

521  IF    (DUR1-DUR2)    7,8,8  J  37 

522  7  IEND<  ID)=DUR2/DT(  IDH-1.  J  38 

523  M3=DUR1/0T(  ID>«-1.  J  39 

524  K1=ID2  J  40 

525  GO   TO  9  J  41 

526  8  IEND(ID)=DUR1/DT(ID>«-1.  J  42 

527  M3=DUR2/DT( ID)*1.  J  43 

528  K1  =  ID1  J  44 

529  9  IF    <  IEND(  ID)-300)    11,11,10  J  45 

530  10         IEND(ID)=300  J  46 

531  11         M2=IEND(K)  J  47 

532  J  =  l  J  48 
C            INTERPOLATE  ONE   HYDROGRAPH    IF   NECESSARY  J  49 

533  TIDH=0.  J  50 

534  T I D  =  DT (  ID )  J  51 

535  DO   15   1=2, M2  J  52 

536  TIDH=T  IDH«-DT(  K)  J  53 

537  12         IF    (TIDH-TID)    15,13,14  J  54 

538  13         J=J+1  J  55 

539  DATA   ( J ) =0C  FS ( I  ,  K )  J  56 

540  T I D=T I D+DT ( ID)  J  57 

541  IF    (J-300)    15,16,16  J  58 

542  14         J  =  J+1  J  59 

543  DATA    (  J  )  =OCFS  (  I  -  1  ,  K  )  +  (  (  T  I  D-T  IDH*-DT<  K  )  )  /  DT  (  K  )  )  *  ( OC  F  S  (  I  ,  K  ) -OC  FS  (  I -1  ,      J  60 
IK  )  )  J  61 

544  T I D=T I D  +  DT (  ID  )  J  62 

545  IF   (J-300)    12,16,16  J  63 

546  15         CONTINUE  J  64 

547  16         IEND(K)=J  J  65 

548  DO   17   1=2, J  J  66 

549  17         OCFS( I , K ) =  DAT  A ( I )  J  67 

550  18         M=IEND(ID»  J  68 

551  DA(ID)=DA( ID1)+0A< ID2)  J  69 

552  R0=0.  J  70 
C            ADC  HYDROGRAPHS  J  7* 

553  DO  2C   1=1, M3  J  72 

554  OCFSI I , ID)=OCFS( I ,ID1 )+OCFS( I , I D2 »  J  73 

555  IF   (OCFSU  ,  I D )  -  PE  AK  ( I D )  )    20,20,19  J  74 

556  19         PE AK( I D)=OCFS( I , ID)  J  75 

557  20         RO=RO+OCFS ( I, ID)  J  76 

558  IF    (PEAM I C)-PEAK(K  1) )    21,22,22  J  77 

559  21         PE AK( I D)=PEAK (Kl )  J  78 

560  22         IF    (M-M3)    25,25,23  J  79 

561  23         M3=M3+1  J  80 

562  DO   24   I=M3,M  J  81 

563  OC  FS  (  I  ,  I D )  =0C  FS  (  I  ,K  1 )  J  82 

564  24         RO=RO+OCFS( I, ID)  J  83 

565  25         ROIN=(RO*DT( ID) )/(DA( ID)* 64 5. 333)  J  84 

566  IF    (NPU)    27,27,26  J  S5 

567  26         WRITE    (7,28)    I D ,  NHD , I D 1 , 1 D2  J  86 

568  27         RETURN  J  87 
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SL0PE1=DATA(6 I 
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CHANNEL  SLCPE. 
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L 

1  a 

c 

Cl     AAA      Ol    ATM  ClAAC 

rLLUD    PLAIN  SLOPE. 

i 

L 

O  A 

20 

595 

ATC/CMAV  CIAI/1A 

DI r=l EMAX  —  ELO  )/ IV. 

L 

0  1 
^  I 

596 

C  (  1  ) =E  L0 

L 

0  0 

2  <: 

59  7 

AA       1       f  —  O  OA 

DO    1  1=2,20 

1 

L 

C  3 

59  8 

1 

l~   1    T  1  —A  IT       11  j-A  T  C 

III) =C I  1  -  I  )  *u  I  r 

L 

0  /. 

c 

p  r  T      ADCA       AMA      ATCALJADAr"       ADOAl/C    A 

SET   AREA   AND    DISCHARGE    ARRAYS   =  U. 

L 

0  K 

599 

r\  r\     o     T        1  OA 

DO    2  1=1,20 

i 

L 

O  A 

AAA 

A  /    T      T  n  1  —  A 

Al  1  f  tU  J  — U  • 

1 

i. 

601 

2 

Q {  I , ID  )=0. 

L 

28 

602 
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WRITE    (6,24)  VS 
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READ  N   VALUES   AND  SEGMENT   BORDER  POINTS. 
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31 
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DO   3  I=1,NSEG 

L 
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605 

SEGNI I )=DATA( J) 

L 

33 

606 

DISK I)=DATA( J+l) 

L 

34 

607 

3 

J  =  J  +  2 

L 

35 

C 

REPAINTING  DATA   ITEMS   ARE   DISTANCES   AND  ELEVATIONS. 

L 

36 

608 

J  J  J=J 

L 

37 

609 

DO  6  I=1,MSEG 

L 

38 

610 

4 

J  =  J+2 

- 

L 

39 

611 

IF    (DATA( J )-D  1ST (  I )  )  4,5,5 

L 

40 

612 

5 

I SG(  I  ) =J+1 

L 

41 

613 

6 

CONTINUE 

L 

42 

C 

COMPUTE   DISCHARGES   AND   END  AREAS 

FOR  EACH 

SEGMENT. 

L 

43 

614 

DO   22   K= 1 , NSE  G 

L 

44 

615 

J=  JJJ 

L 

45 

61  6 

J  J  Jl = J  J J+ 1 

L 

46 

617 

IF    (SEGN(K))  7,7,8 

L 

47 

618 

7 

SL0PE=SL0PE1 

L 

48 

61  9 

S  E  GN ( K ) =-S  EGN ( K  ) 

L 

49 

620 

GO   TO  9 

L 

50 

621 

8 

SL  CPE=SL0PE2 

L 

51 

622 

9 

SLPN=1. 48 6* SLOPE**. 5 

L 

52 

C 

COMPUTE   AREA  AND   DISCHARGE   FOR  SEGMENT. 

L 

53 

62  3 

DO   21  1=2,20 

L 

54 

624 

AA=0. 

L 

55 

62  5 

P  =  0. 

L 

56 

626 

J= JJJ-1 

L 

57 

627 

DEP2=0. 

L 

58 

628 

10 

J  =  J  +  2 

L 

59 

629 

IF    (J-ISG(K))  12,12,11 

L 

60 

630 

1  1 

IF    (AA-.OOl)  21*21,20 

L 

61 

631 

12 

IF    (DATA! J)-Ct I ) 1  13,10,10 

L 

62 

632 

1  3 

DEP1  =  C(  I  ) -D AT  A ( J ) 

L 

63 

633 

IF    (J-JJJ1)  16,16,14 

L 

64 

634 

14 

XL=DATA( J-l )-DATA( J-3 ) 

L 

65 

63  5 

D  E  P3  =  A  B  S ( D  A  T A ( J-2 )-CATA( J ) ) 

L 

66 

636 

XL=XL*DEP1/DEP3 

L 

67 

637 

15 

AA=AA+XL*( DEP1+DEP2 )/2. 

L 

68 

638 

P=P+SQPT((DEP1-DEP2)**2+XL**2) 

L 

69 

639 

16 

DEP2=DEP1 

L 

70 

640 

J  =  J+2 

L 

71 

641 

IF    (J-ISG(K))  17,17,20 

L 

72 

642 

1  7 

IF    ( DATAI J )-C ( I » )  18,18,19 

L 

73 

64  3 

1  8 

DEP1=C( I )-OATA( J) 

L 

74 

644 

XL=DATA( J-l  ) -DATAI J-3) 

L 

75 

645 

GO   TO  15 

L 

76 

646 

19 

D  E  P 1 =0  . 

L 

77 

64  7 

XL=DATA(J-1)-DATA( J-3) 

L 

78 

64  8 

DEP3=ABS(DATA( J-2)-DATA(J) ) 

L 

79 

64  9 

XL=XL*DEP2/DEP3 

L 

80 

650 

AA=AA+XL*( DEP1+DEP2 )/2. 

L 

81 

65  1 

P=P+SQRT((DEP1-DEP2)**2+XL**2) 

L 

82 

652 

DEP2=0  . 

L 

83 

653 

GO   TO  10 

L 

84 

654 

20 

R  =  AA/P 

L 

85 

655 

SGN=SEGN( K )-.0025*R 

L 

86 

C 

ADD  DISCHARGES  AND  AREAS  FOR  ALL 

SE  GMENTS 

TO  OBTAI N  TOTALS  FOR 

L 

87 

C 

VALLEY  SECTION. 

L 

88 

656 

Q(  I,ID)=3<  I,ID)+AA*R**.66667*SLPN/SGN 

L 

89 

657 

A(  I ,  ID  )  =A(  I  ,  ID)  4-AA 

L 

90 

658 

21 

CONTINUE 

L 

91 

659 

JJ J=J-3 

L 

92 

67 


66  0 

c  c 

r  n mt t  mi ic 

UUI\  1  I  IN  U  u 

661 

00   23  1=1,20 

662 

DE  EP ( I  ,  ID ) =C( I )-EL0 

663 

WRITE   (6,25)  C(I),A(I, 

ID)  fQ( I , ID) 

664 

23 

CONTINUE 

665 

r 

RETURN 

666 

\* 

24 

FORMAT ( 1H0 ,T42 , 'RATING 

CURVE   VALLEY   SECTION' ,F5. 1/T46, 'WATER' ,  T56  , 

l'FLOW'  ,T66,'FL0W'/T45, 

•  SURFACE' ,T56,'AREA'  , T66 , • R ATE • /T46 , • EL E V • , T 

256,' SQ   FT • ,T66, 'CFS' ) 

667 

25 

FORMAT    (40X,F  10.2.2F10 

.  1) 

668 

END 

669 

SUBROUTINE   ST  T 

670 


671 
672 
673 
674 
675 
676 
677 
678 
679 
680 
681 
682 
683 

684 


685 


THIS  SUBROUTINE   STORES   A  DEPTH  -  FLOW  -  TRAVEL  TIME  TABLE. 

COMMON  CFS (300) ,OCFS( 300,6) , I  END ( 6  )  , DAT A ( 3 10 > , DA ( 6 ) , DP ( 20 ) , NPU ,NHD 
It SCFSl20)t CI20) ,A(20t6)  ,  Q ( 20 , 6 ) , R A  I N ( 200 ) , DEEP ( 20 , 6 ) , NER, MAXNO, NCO 
2MM,ICC,NC0DE,DIST(6),SEGN(6),CTBLE(50,11) , I TBL E ( 5 0 , 2 ) , Z ALF A ( 20 ) , DT 
31 6 ) t TIME  f PEAK( 6 ),R0IN,ISG(6) 

I D  =DAT A (  1  ) 

R  E  ACH=  DATA  (  2) 

XL  =DAT A ( 3 ) 

SL0PE=DATA(4) 

D I ST( ID)=SLOPE*XL 

J  =  5 

DO    1  1=1,19 

DP ( I )  =  D AT  A ( J ) 

SCFS( I ) =DATA( J+l  ) 

C (  I ) =DAT A ( J  +  2  ) 

J  =  J  +  3 

RETURN 

END 

SUBROUTINE  CMPTT 

THIS  SUBROUTINE  COMPUTES   THE  TRAVEL    TIME  AT  GIVEN 
DISCHARGE  RATES 

COMMON  CFS (300) ,OCFS( 300 , 6 ) , I  END ( 6)  , DAT A ( 3 1 0 ) , DA ( 6 ) , DP ( 20 ) , NPU ,NHD 
It  SCFSf  20),C(20)tA(20f6),Q(20,6) , RAI N ( 200 ) , DEEP ( 20 , 6 ) , NER, MAXNO, NCO 
2MM,ICC,NC0DE,DIST(6),SEGN(6) , CT BLE ( 50 , 1 1 ) , ITBLE(50,2) , ZALFA ( 20 ) , DT 
3( 6),TIME,PEAK(6),R0IN,ISG<6) 


686 

I D=DAT  A ( 1 ) 

687 

RE  ACH=DATA (2) 

688 

N0VS=DATA( 3) 

689 

X  L=DAT A ( 4 ) 

690 

SLOPE=DAT A ( 5 ) 

691 

D I ST( ID)=SLOPE*XL 

692 

XLD36=XL/3600. 

C 

ZERO  ARRAYS 

693 

DO    1  J=l,20 

694 

DATA  (J)=0. 

695 

1 

CFS( J)=0. 

696 

ID1=1 

C 

FIND  RATING  CURVE  WITH  SMALLEST   MAXIMUM  FLOW  RATE 

697 

2 

QM  IN=QI 20, ID1 ) 

698 

MIN-ID1 

699 

GO    TO  4 

700 

3 

I  Dl=  ID  1-t-l 

701 

IF    (QMIN-Q(20, ID1) )  4,4,2 

702 

4 

IF    (ID1-N0VS)  3,5,5 

703 

5 

1=1 

L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 

M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 
M 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 


68 


C 

SET   SCf-S    ARRAY   EQUAL   TO  Q  ARRAY   OF   LOWEST   RATING  CURVF 

N 

28 

704 

00   6  J=2,2C 

N 

29 

705 

SCFS(  I  )=Q( J,MIN) 

N 

30 

706 

6 

1  =  1+1 

N 

31 

C 

COMPUT    END  AREA  AND  DEPTH 

N 

32 

707 

DO   9  ID1=1,N0VS 

N 

33 

708 

DO  9  J =1,19 

N 

34 

709 

DO   7  1=2,20 

N 

35 

710 

IF    ( Q (  I  1 1 D  1 )  — SCFS ( J  )  )  7,17,8 

N 

36 

71 1 

7 

CONTINUE 

N 

37 

712 

17 

DA  TA   ( J  )  =  A (  1,10  1) +  D AT  A ( J ) 

N 

38 

713 

CFS(J)=DEEP(I,ID1)*CFS(J) 

N 

39 

714 

GO   TO  9 

N 

40 

715 

8 

XY  =  ( SCFS( J  )-Q(  1-1, ID1 ) )/( C( I , ID1  )-Q( 1-1, ID  1)1 

N 

41 

716 

DATA   (  J)=A(  I-  1,  ID1)+XY*(  A(  I  ,  IDi)-A(  1-1,  ID1  > )+DATA(  J) 

N 

42 

717 

CFS(J)=DEEP(I-l,lDl)+XY*(DEEP(I, I D  1 ) - D  E  E  P (  I- 1 ,  I  D 1 )  >  +C  F  S  (  J  ) 

N 

43 

71  8 

9 

CONTINUE 

N 

44 

719 

XNCVS=NOVS 

N 

45 

720 

WRITE   (6,13)  REACH 

N 

46 

C 

COMPUTE    TRAVEL  TIME 

N 

47 

721 

DO    10  1=1,19 

N 

48 

722 

A V  AR  EA  =  D AT  A ( I  )/XNOVS 

N 

49 

72  3 

DP  (  I  )  =  CFS(  I  )/XNOVS 

N 

50 

724 

S=AVAREA*XLD36 

N 

51 

72  5 

C (  I )  =S/SCFS(  I  ) 

N 

52 

726 

WRITE    (6,14)    DP ( I ) , SC  FS ( I ) , C ( I  ) 

N 

53 

727 

10 

CONTINUE 

N 

54 

c 

PUNCH  CODE 

N 

55 

728 

IF    (NPU)  12,12,11 

N 

56 

729 

11 

WRITE   (7,15)    ID, REACH, XL, SLCPE 

N 

57 

73  0 

WRITE    (7,16)    ( DP(  I ) ,SCFS( I »  ,C( I  )  ,  1=  1, 19) 

N 

58 

731 

12 

RE  TURN 

N 

59 

c 

N 

60 

732 

1  3 

F □ RM AT ( 1H0 ,T46,' TRAVEL   TIME    TABLEVT54, 'REACH' ,F5. 1//T46, 'WATER'  ,T 

N 

61 

156 , 'FLOW'  , T65 ,' TRAVEL ' /T46,  'DEPTH'  ,  T 56 , ' R A TF •  , T66 , • T I  ME • / T46 » ' FE E T 

N 

62 

2'  , T56,  'CPS'  ,T66» ' HR S ' ) 

N 

63 

733 

14 

FORMAT  (40X,F10.2,F10.0,F10.4) 

N 

64 

734 

15 

FORMAT  (  'STORE    TRAVEL    TI      •  ,  T21  ,  •  I  D=  •  ,  1 1  ,  T2  9,  •  REACH  NO  =  '  ,  F  5.  1  ,  T44  , 

N 

65 

1 « L ENGTH= • , F9.0 ,  •   F T • / T 2  1 ,  '  S LOP E= • , F 8. 6 , • F T /FT •  /T2 

N 

66 

21 ,  'DEPTH! FT)'  ,T35,'FL0W(CFS)'  ,  T49 , • T I MF ( HRS ) *  1 

N 

67 

73  5 

16 

FORMAT   (T2i,F7.2tF15.0,F15.  3) 

N 

68 

736 

END 

N 

6  9- 

737 

SUBROUTINE  ROUTE 

0 

1 

C 

THIS   SUBROUTINE   ROUTES   A  HYDROGRAPH   THROUGH   A   REACH  WITH  THE 

0 

2 

C 

NEW  VSC   METHOD  OF  FLOOD  ROUTING.    THIS  METHOD  ACCOUNTS   FOR  THE 

0 

3 

C 

VARIATION    IN   WATER   SURFACE  SLOPE. 

0 

4 

738 

COMMON   CFS  (300)  ,OCFS( 300,6),IEND(6),DATA(310),DA(6),DP(20) ,NPU,NHD 

0 

5 

1,SCFS(20),C(20),A(20»6),Q(20,6), R A  I N ( 200 ) , DEEP ( 20 , 6 ) , NER , MA XNO , NCO 

0 

6 

2MM  ,1 CC,NCODE,D  I  ST (6 ) , SEGN (6) ,CTBLE( 50, 11 ) , ITBLE(50 ,2) , Z ALF A ( 20 ) , DT 

0 

7 

3(6),TIME,PEAK(6),R0IN,ISG(6) 

0 

8 

739 

I D  =  DAT  A ( 1 ) 

0 

9 

740 

NHD=DATA( 2 ) 

0 

10 

741 

I DH  =  DA  T A ( 3 ) 

0 

1  1 

742 

DT ( I D ) =DAT  A ( 4 ) 

0 

12 

743 

DA ( ID)=DA( IDH) 

0 

13 

744 

M=  I  END (  IDH) 

0 

14 

C 

IF    ID  AND    IDH   ARE   EQUAL,   ADO   1    TO  IDH 

0 

15 

745 

IF    (ID-IDH)  3,1,3 

0 

16 

746 

1 

I0F=IDH+1 

0 

17 

69 


747 

748  2 

749 

750 

751  3 

752 

753 

754 

755 

756 

757 

758 

759 

760 

C 
C 

761 

762  4 

763 

764 

765 

766 

767  5 
76  8 
769 
770 

771  6 
772 


773 
774 
775 
776 
777 
778 
779 
780 
781 
782 
783 
784 
785 
786 


787 
788 
789 
790 
791 
792 
793 
794 

795 
796 
797 

798 
799 
800 


9 

10 


1  1 


12 
13 


14 

C 

1  5 

16 

C 

49 


EQUAL  TO  TIME   INCREMENT   OF  INFLOW 


00  2    1  =  1, M 

OCFS( I , IDH) =OCFS( I , IDH-1) 
DT( I  DH ) =  DT ( IDH-1 ) 
P  E  AK(  IDH)  =  PEAK< IDH-1) 
NE  PRT=0 
PE  AK( I D)  =  l  . 
R0=0. 
N=  19 

OCFSd  ,  ID)=0. 

s=o. 

T1=C( 1 ) 
J=l 

GUES=1 . 
CFS(  1)  =0. 

IF   ROUTING   INTERVAL   IS  NOT 
HYDROGRAPH,  INTERPOLATE 
IF    ( DT(  IO)-DT(  IDH)  )  8,15,4 
T I D=  DT ( ID) 
T I DH=0. 
DO   7   1=2, M 
T I DH=T IDH+DT ( IDH) 
IF    (TID-TIDH)  6,5,7 
J  =  J+l 

CFS(  J  )  =OCFS(  I  ,  IDH) 
TIC=TID+DT( ID) 
GO  TO  7 
J  =  J+l 

CF  S( J)=OCFS( 1-1, IDH)  +  (  ( TIO-TIDH+DTI IDH) )/DT( IDH) ) *< OCFS ( I , IDH)-OCF 
1S(  1-1,  IDH)  > 
TIC=TID+DTUD) 
CONTINUE 
GO   TO  13 
TI DH=0. 
T I D=DT{ ID) 
DO   12   1=2, M 
T I DH=T IDH+DT (IDH) 
IF    (TIDH-TID)  12,10,11 
J  =  J+l 

CF  S(  J)  =OCFS(  I  ,  IDH) 

TIC=TID+DT(ID) 

IF    (J-300)  12,13,13 

J=  J+l 

CF  S( J) =OCFS( 1-1 » IDH)* ( ( TI D— T IDH+DT (  I  DH ) > /DT ( I DH ) ) * ( OCF S ( I , I DHJ-OCF 
1S(  I- 1, IDH)  ) 
T I C=T I D+DT ( ID ) 
IF    (J-300)  9,13,13 
CONTINUE 
IENDl IDH)=J 
DT (  IDH)=UT( ID) 
M=  J 

DO    14   1=2, ^ 

OCFS( I , IDH)=CFS(  I  ) 

IF    INFLOW    IS    ZERO ,    SO   IS  OUTFLOW 
DO    15  L=2,M 

IF    (OCFSIL, IDH) )  16,16,49 

OCFS(L  ,  ID) =0. 

ROUTE 

DATA  ( L-l ) =0. 
DO  42  I=L,300 
IF    (I-M)  18,18,17 


0 

18 

0 

19 

0 

20 

0 

21 

0 

22 

0 

23 

0 

24 

0 

25 

0 

26 

0 

27 

0 

28 

0 

29 

0 

30 

0 

31 

0 

32 

0 

33 

0 

34 

0 

35 

0 

36 

0 

37 

0 

38 

0 

39 

0 

40 

0 

41 

0 

42 

0 

43 

0 

44 

0 

45 

0 

46 

0 

47 

0 

48 

0 

49 

0 

50 

0 

51 

0 

52 

0 

53 

0 

54 

0 

55 

0 

56 

0 

57 

□ 

58 

0 

59 

0 

60 

0 

61 

0 

62 

0 

63 

0 

64 

0 

65 

0 

66 

0 

67 

0 

68 

0 

69 

0 

70 

0 

71 

0 

72 

0 

73 

0 

74 

0 

75 

□ 

76 

0 

77 

70 


801 

1  7 

0CFS( I  ,  IDH>=0CFS(  I- 1, IDH)*.9 

0 

78 

802 

18 

AVlN=(OCFS( It  IDH)*OCFS(  I-l,  IDH)  )/2. 

0 

79 

803 

SI A=S+AVIN 

0 

80 

804 

J  =  l 

0 

81 

C 

DETERMINE    DEPTH  AND   TRAVEL   TIME  OF 

INFLOW 

0 

82 

805 

IF    { OCF S ( I  ,IDH)— SCFSI  I)  )  19,23,20 

0 

83 

806 

19 

U I 2=(0CFS<  I,IDH)/SCFS(D)  *DP(  1  ) 

0 

84 

807 

T I 2=C(  1 ) 

0 

85 

808 

GO  TO  25 

0 

86 

809 

20 

DO   21  J=2,N 

0 

87 

810 

IF    (  0  C  F  S  (  I  ,IDH|  —  SCFS(J)  )  24,23,21 

0 

88 

811 

21 

CONTINUE 

0 

89 

812 

IF    (NERRT)  22,22,36 

0 

9C 

813 

22 

WRITE  (6,461 

0 

91 

814 

NE  RR  T= 1 

0 

92 

815 

GO   TO  36 

0 

93 

816 

23 

DI  2-DP  (J) 

0 

94 

817 

T I 2  =  C<  J  ) 

0 

95 

818 

GO  TO  25 

0 

96 

819 

24 

RATIO=(OCFSU,IDH)-SCFS(J-l))/<  SCFS  <  J  )- 

SCFS(J-l) ) 

0 

97 

820 

DI2=DP(J-l )+RATIO*(DP( J)-DP(J-l) > 

0 

98 

821 

TI 2=C< J-l) «-RATIO*(C (J)-C( J- 1 ) ) 

0 

99 

822 

2  5 

DO   35  IT=1,1C 

0 

ion 

823 

J  =  l 

0 

101 

C 

DETERMINE    DEPTH  AND   TRAVEL   TIME  OF 

OUTFLOW 

0 

102 

824 

IF    ( GUE  S-S  CFS ( I ) )  26,29,27 

0 

103 

825 

26 

D0  2=(GUES/SCFS<  1) )*DP( 1 ) 

0 

104 

826 

T02=C( 1 > 

0 

105 

827 

GO  TO  31 

0 

106 

828 

27 

DO   28  J=2,N 

0 

107 

829 

IF    ( GUES-SCFS ( J ) )  30,29,28 

0 

108 

830 

28 

CONTINUE 

0 

109 

831 

J  =  N 

0 

110 

8  52 

29 

D02=DP ( J ) 

0 

1  1  1 

833 

T02=C( J) 

0 

112 

834 

GO  TO  31 

0 

113 

835 

30 

RATIO=(GUES-SCFS(  J- 1 )  )  /  (  SCFS  {  J  )-SCF  S  U- 

1)  ) 

0 

114 

836 

D02=DP< J-l )+RATIO*(DP( J)-DP(J-l) ) 

0 

115 

837 

TO  2  =  C ( J-l ) ♦RATIO*(C ( J )-C( J-  1 ) ) 

0 

116 

C 

FIND  WATER    SURFACE  SLOPE 

0 

117 

838 

31 

DD  D=D I S  T (  IC)/(DIST( I D ) ♦ C I 2-D02 ) 

0 

118 

839 

IF    (DDD-.01)  32,32,33 

0 

119 

840 

32 

GUES=0CFS(  I-l  ,  I DH  » 

0 

120 

841 

GO  TO  35 

0 

121 

842 

33 

T2=.5*( TI2+TU2) 

0 

122 

843 

T2=T2*SQRT (ODD) 

0 

123 

844 

T=T1+T2 

0 

124 

C 

COMPUTE   ROUTING  COEFFICIENT 

0 

125 

845 

C0EF=(2.*DT( I  0 ) ) / 1 T  *DT { I  0 ) > 

0 

126 

84  6 

02=C0EF*S I  A 

0 

127 

84  7 

TR  Yl =GUES 

0 

128 

84  8 

RATI  0=02/1 GUES+.1E-20) 

0 

129 

849 

D  I  FF=ABS(  1  .-RATIO 

0 

130 

C 

TEST  FOR  CONVERGENCE 

0 

131 

850 

IF    (DIFF-.001 )  37,37,34 

o 

132 

851 

34 

GUES=02 

0 

133 

852 

35 

CONT  INUE 

0 

134 

853 

OC  FS( I , ID) =DATA( I-l )*Sl A 

0 

135 

854 

DATA   (  I  )  =D AT A ( I-  1 ) 

0 

136 

855 

WRITE    ( 6,47  >    I  , OC  F  S ( I ,  ID) 

0 

137 

71 


856  GO   TO   38  0 

857  36  OCFS( I , ID) =DATA( 1-1 )*S  I  A  0 

858  DATA   (  I ) =D A TA ( 1-1)  0 

859  GO   TO   38  0 

860  37  OC  FS( I f ID ) =02  0 

861  DATA  ( I )=COEF  0 
C  COMPUTE   NEW   STORAGE  0 

862  38  S=SIA-OCFS( I, ID)  0 

863  T1=T2  0 

864  RO=RO+QCFS(  I »  ID)  0 

865  IF    ( OCFS< I  ,  ID l-OCFS ( 1-1  ,  ID) )    39,40,40  0 

866  39  IF    (OCFS(  I  ,  ID  )-l.  )    43,43,42  0 

867  40  IF    ( OCFSI  I  ,  ID  )-PE AK ( I D)  )    42,42,41  0 

868  41  PEAK( ID)=OCFS ( I  ,ID)  0 

869  42  CONTINUE  0 

870  1  =300  0 

871  43  IEND(ID)=I  0 

872  ROIN=(RO*DT( ID)  )/(DA(  ID)* 64 5. 333)  0 
C  PUNCH  CODE  0 

873  IF    (NPU)    45,45,44  0 

874  44  WRITE   (7,48)    I D , NHD , I DH , DT ( I D )  0 

875  45  RETURN  0 
C  0 

876  46  FORMAT! 1H0 , 'TRAVEL   TIME  TABLE    EXCEEDED')  0 

877  47  FORM  AT (T10»* PROBLEM  FAILED  TO  CONVERGE   AFT  ER 10    ITERATIONS.   CONVERG  0 

1ENCE  WAS  FORCED. '/T20 , 'OUTFLOW  NUMBER  =   ',14, 'RATE   =',F10.2)  0 

878  48  F  0  RM  AT (        'ROUTE', T21, '10  =  ', I1»T29,'HYD  N0=*, 13, T45, 'INFLOW    ID=',I  0 

11 ,T65,'DT=',F8.6,'HRS' )  0 

879  END  0 

880  SUBROUTINE  RESVO  P 
C  THIS  SUBROUTINE  ROUTES  A  HY DROGRAPH  THROUGH  A  RESERVOIR  WITH  THE  P 
C  STORAGE-INDICATION  METHOD.  P 

881  COMMON  CFS ( 300)  ,OCFS( 300,6)  , IEND(6  )  ,DATA( 310) ,DA( 6) , DP( 20) , NPU, NHD  P 
1,SCFS(20),C(20),A(20,6),Q(20,6) ,RAIN(200) , DEEP (20, 6) , NER , MA XNO , NCO  P 
2MM,ICC,NC0DE,DIST(6),SEGN(6) , CT BLE ( 50 , 1 1 ) , ITBLE(50,2) , ZALF A ( 20 ) , DT  P 
3( 6),TIME,PEAK(6) ,R0IN,ISG(6)  P 

882  I D  =  DAT  A ( 1 )  P 

883  NH  D=  DA  T  A ( 2  )  P 

884  I DH=DA  T A (  3  )  P 

885  N  E  RE  S=0  P 

886  DT ( ID)=DT (  IDH )  P 

887  R0=0.  P 

888  DA ( ID)  =  DA(  IDH  )  P 

889  PEAK(ID)=1.  P 

890  J=l  P 

891  1=4  P 
C  REMAINING  DATA   ARE   FLOW  AND  STORAGE  VALUES  P 

892  SC  FS ( J  >  =DAT A (  I )  P 

893  STCRE1=DATA(I*1)*12.1  P 

894  ST0RE=ST0RE1  P 
C  COMPUTE   STORAGE  COEFFICIENT  ARRAY  C  P 

895  1  C ( J)  =  ( SCFS( J >/2. )  +  ( STORE/DT( ID) )  P 

896  1=1+2  P 

897  J=J+1  P 

898  IF    (J-20)    2,2,3  P 

899  2  SC  FS ( J ) =DA  T A (  I  )  P 

900  STORE=DATA ( I  +  l)*12.  1  P 

901  IF    (SCFS(J)-.OOl)   3,3,1  P 

902  3  N=J-1  P 
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91  5 

916 

8 

c 

91  7 

918 

91  9 

9 

920 

921 

1  0 

c 

922 

1 1 

92  3 

92  4 

925 

12 

92  6 

1  3 

927 

1  4 

92  8 

1  5 

92  9 

930 

16 

931 

C 

932 

93  3 

I  7 

93  A 

93  5 

'  -J '  j 

1  8 

C 

937 

19 

938 

20 

939 

2  1 

940 

941 

C 

C 

942 

94  3 

944 

94  5 

94  6 

947 

□CFS(l,IO>=0. 
S=  STDR  E1/DT(IC) 
C  ROUTE 

00   15  1=2,150 

IF    (I-  IEND( IDH)  )    5, 5,4 

4  OC  FS (  I  »  I  OH  )  =0 . 0 

5  AVIN=(OCFS(  I,  IDH) +OCFS (  1-1  ,  IDH)  ) /2. 
51  A=S+-AVIN 

C  DETERMINE    PROPER  C 

DO   6  J=1»N 

IF   (SIA-C(J))  10,9,6 

6  CONTINUE 

IF    (NERES)  7,7,8 

7  WRITE  (6,19) 
N  E  RE  S=  1 

RESC=SCFS(N)/C(N) 
COMPUT  OUTFLOW 
OCFS( I , ID) =RESC*S  IA 
GO  TO  11 

OC  FS( I , ID) =SCFS ( J ) 

GO  TO   11  , 
OCFS(  I , ID)=SCFS<  J-l )>( (SIA-C(J-l) )/ (C(J)-C (J-l) > )* (SCFS( Jl-SCFSl  J- 

11)  ) 

DETERMINE   NEW  STORAGE 
S=SI  A-OCFS(  I,  ID) 
RO=RO«-OCFS  (  I  ,  ID) 

IF    (OCFS(  I ,  ID )-OCFS ( 1-1,10) )  12,13,13 
IF    (OCFS(  I  ,  ID  )-l.  )  16,16,15 
IF    (OCFSI I , ID )-PEAK ( I D) )  15,15,14 
PEAK(ID)=OCFS(I ,ID) 
CONTINUE 
I  =  150 

I  F  ND  (  I  D  )  =  I 

ROIN=RO*DT( ID)/(DA( ID)*64  5.333) 
PUNCH  CODE 
IF    (NPU)  18,18,17 
WRITE    (7,20)    I  D  ,NHD , I DH 
I  I =2*N+3 

WRITE    (7,21)    (DATA( I )  ,  1  =  5, I  I ) 
RE  TURN 

FORMAT    ( 1H0 , 3 3HST0R AGE-D I SCHARGE   TABLE  EXCEEDED.) 

F  0  PM AT (        'ROUTE   RESERVOIR* ,T21»* ID=*  tIl,T29,  *HYD  N0= ',13  »T42 ,  ' INF 
1L0W   ID=',U  /T21,  'OUTFLOW(CFS)  •  ,T37,  •  STOR 

2AGE(  AC   FT  )  •  ) 
FORMAT    ( T21  ,F10.1,F 13. 1) 

END 

SUBROUTINE  ERROR 

THIS  SUBROUTINE   DETERMINES  THE   ERROR  STANDARD   DEVIATION  AND  THF 
PEAK  FLOW   ERROR  FOR   2  HYDROGRAPHS 

COMMON   CFS (300)  ,OCFS( 300,6)  ,  I  END (6)  , DAT A ( 310) , DA ( 6 ) , DP ( 20 ) , NPU ,NHD 
1,SCFS( 20)  ,C(20)  ,A(20,6)  ,Q(20,6)  , RAIN (200) , DEE P ( 20 , 6 ) , NFR , MA XNO , NC 0 
2MM,ICC,NCODE,DIST(6), SEGN(6) , CT BLE ( 50 , 1 1 ) , ITBLE(50,2) ,ZALFA(20) ,DT 
3( 6  ) , TIME.PE AK(6) ,ROIN, ISG(6) 
I D  1=DATA(  1 ) 
I D2=0ATA( 2 ) 
SSE=0. 

WRITE  (6,15) 
J=l 
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C 

IF   TIME   INCREMENTS  NOT  EQUAL,  INTERPOLATE 

Q 

13! 

94  8 

IF   (DT ( IDl l-OT( ID2I )  1,8,2 

Q 

14 

949 

1 

L=ID1 

Q 

15!) 

95C 

K=  ID2 

Q 

16 

951 

GO  TO  3 

0 

17 

952 

2 

L=  ID2 

Q 

18 

953 

K=  101 

Q 

19| 

954 

3 

M=IEND(L> 

Q 

201 

955 

T I D=DT ( K ) 

Q 

21 

95  6 

TIDH=0. 

Q 

22  1 

957 

00  6  1=2, M 

Q 

23 

95  8 

T  I  DH=T IDH+DTt  L) 

Q 

24 

959 

IF   (TID-TIOH)  5,4,6 

Q 

25 

960 

4 

J=  J+l 

Q 

26 

961 

£FS(  J)=OCFS(I  ,1_L 
T I D=T  I  D+DT ( K)  X 

Q 

27 

96  2 

Q 

28 

963 

GO   TO  6 

Q 

29 

964 

5 

J  =  J+l 

Q 

30 

96  5 

CFS(JI=0CFS(I-1,L)+((TID-TIDH*DT(L) )/CT(L) )*(OCFS( I, 

L)-0CFS(I-1,L) 

Q 

31 

1) 

Q 

32 

966 

TIOTID  +  DT (K> 

Q 

33 

96  7 

6 

CONT  INUE 

Q 

34 

968 

I END(L ) =J 

Q 

35 

969 

DT(L)=DT(K  ) 

Q 

36 

970 

DO   7   1=2, J 

Q 

37 

971 

7 

OC  FS( I ,  L ) =CFS  (  I  ) 

Q 

38 

972 

8 

IF    (IEND( ID1)-IEND( ID2) )  9,9,10 

Q 

39 

973 

9 

M=  I  END ( IDl ) 

Q 

40 

974 

GO  TO  11 

Q 

41 

97  5 

10 

M=  I  END ( ID2 ) 

Q 

42 

976 

11 

T  2  =T I  ME 

Q 

43 

C 

DETERMINE  ERROR 

Q 

44 

977 

DO   12    1  =  1,  M 

Q 

45 

978 

ERR=OCFS( I ,ID1)-0CFS( 1,102) 

Q 

46 

979 

WRITE    (6,16)    T2,0CFS( I, IDl) ,OCFS(I, 

ID2) ,ERR 

Q 

47 

980 

T2  =  T2+DT(  IDl) 

Q 

48 

C 

SUM  OF    SQUARES  OF  ERROR 

Q 

49 

98  1 

12 

SSE=SSE*ERR*ERR 

Q 

50 

982 

XM=M 

Q 

51 

C 

ERROR  VARIANCE 

Q 

52 

983 

EVAR=SSE/XM 

Q 

53 

C 

ERROR  STANDARD  DEVIATION 

Q 

54 

984 

ESDEV=SQRT(EVAR) 

Q 

55 

985 

WRITE   (6,17)  ESDEV 

Q 

56 

C 

PERCENT   ERROR   FOR  PEAK  FLOWS 

Q 

57 

986 

ERPK=ABS(PEAK(ID1)-PEAK(ID2) » 

Q 

58 

987 

PCTER=(ERPK/PEAK(ID1) )*100. 

Q 

59 

988 

WRITE   (6,18)  PCTER 

Q 

60 

C 

PUNCH  CODE 

Q 

61 

989 

IF    (NPUI  14,14,13 

Q 

62 

990 

13 

WRITE    (7,19)  101,102 

Q 

63 

99  1 

14 

C 

RETURN 

Q 
Q 

64 
65 

992 

1  5 

FORMAT( 1H0 ,T3  3, 'TIME'  ,T55 , 'FLOW  1', 
l'HRS*  ,T57,  'CPS'  ,T78, 'CPS' t  T97, 'CPS' 

T76, 'FLOW  2'  ,T95, 
) 

•ERROR' /T34, 

Q 
Q 

66 
67 

99  3 

16 

FORMAT  (20X,F20.3,3F20.0) 

Q 

68 

994 

17 

FORMAT! IHOiTIO, 'ERROR   STANDARD  DEVIATION  =  ',F10.3> 

Q 

69 

995 

18 

F0RMAT(T10,'PEAK  DISCHARGE    ERROR  = 

',F7.2,'  PERCENT' 

///) 

Q 

70 

996 

19 

FORMAT (        'ERROR   ANALYSIS' , T21, ' ID 

1=' , 1 1, T29, • ID   I  I 

=',I1> 

0 

71 

997 

ENC 

Q 

72- 
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998 
999 


1000 

1001 

1002 

10C  3 

1004 

C 

1005 

1006 

1007 

C 

1008 

1009 

1 

1010 

2 

C 

ion 

3 

1012 

4 

1013 

SUBROUTINE  SEDT 

THIS  SUBROUTINE  COMPUTES  THE  SEDIMENT  YIELD  FOR  A  FLOOD 

COMMON  CFS (300) ,OCFS( 300,6)  , I  END ( 6 > , D ATA ( 3 10 ) , DA ( 6 ) , DP ( 20 » , NPU , NHD 
I, SCFS( 20),C(20) ,A(20,6) ,Q ( 20 , 6 ) , RAI N< 200 > , DEEP ( 20 , 6 ) ,NER , MAXNO, NCO 
2MM,ICC,NC0DE,DIST(6),SEGN(6),CTBLE(50,11> , I TBL E ( 50 , 2 ) , Z ALF A ( 20 ) , DT 
31 6),TIME,PEAK(6) ,ROIN, ISG(6) 

I 0=DATA( 1) 

SO  IL  =  DATA( 21 

CROP  =  DATA<  3) 

CP=DATA(4) 

SL=DATA( 5) 

COMPUTE   SEDIMENT  YIELD 

X=ROIN*DA( ID)*53.33  3*PEAK(IDI 

SED=95.*X**.56*S0I L*CROP*CP*SL 

WRITE   (6,3)  SED 

PUNCH  CODE 

IF   (NPU)  2,2,1 

WRITE   (7,4)    ID, SOIL, CROP, CP,SL 
RE  TURN 


FORMAT    (10X,    'SEDIMENT   YIELD  =   '.  F10.1,  « 
F  ORM AT (        'SEDIMENT  YIELD'  ,T2l ,' ID=' , II ,T29 
1='  ,F5.3,T57,'CP  =  ' .F5.3.T70, 'LS=' ,F5.3) 

END 


=   ',  F10.1,    '  TONS') 

l,iIO=' , II ,T29,' SOIL=' ,F5.3,T42, 'CROP 

c  _  •     c  c  at 
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24- 


//iDATA 
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ZALFA  =    1234567890  *. ,- 


COMMAND  TABLE 


START  1  2 

STORE  HYO  2310 

RECALL   HYD  3310 

COMPUTE    HYD  4310 

PRINT  HYD  5  2 

PUNCH  HYD  6  1 

PLOT  HYD  7  2 

ADD  HYD  8  4 

STORE  RATING  CURVE  9100 

COMPUTE    RATING  CURVE  10310 

STORE  TRAVEL  TIME  11100 

COMPUTE    TRAVEL  TIME  12  5 

ROUTE  13  4 

ROUTE  RESERVOIR  14100 

ERROR  ANALYSIS  15  2 

SEDIMENT  YIELD  16  5 

FINISH  17  0 
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